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SECTION  I. 

d  - 
d  = 

Ni- 

\/and  £  = 

R  = 


diameter  of  droplet 

mean  diamc-ter  of  spray 

number  of  droplets  diameter  di 

total  volume  and  surface  of  N  droplets 

volume  fraction  of  spray  having  diameter^  d 


Tl  = 

Vb  = 

k  = 

Sivj.Klyel'c  = 


V  = 


size  constant  in  Rosin  Hammier  equation 

distribution  constant  in  Rosin  Rommler 
equation 

volume  rate  of  burning  of  propellant 

proportionality  constant 

surface,  number  etc.,  per  unit  volume 
of  liquid 

volume  rate  of  input  of  propellant 
volume  fraction 


P- 

t  = 
M  - 


index  in  the  equation  relating  rate 
of  burning  to  diameter  of  drop 

time  (seconds) 

rate  of  change  or  flow  of  mass 


TTl  —  H  if  if  [f  it  ,  , 

Tr  “  per  unit  area. 

h  =  thermal  c  c  n  d  uc  t i v i t y 

Up  =  specific  heat  at  constant  pressure 
0  =  non-dimensional  temperature  defined  in  text 
T  -  temptrnturc  in  °A 
L  =  heat  of  vaporisation 
D  =  diffusion  conet“nt 
V  =  kinematic  viscosity 
=  flow  velocity 
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SECTION  II 


r  = 

density 

V  = 

surface  tension 

representative  lengths 

/  - 

viscosity 

Q  = 

volume  flow  rate 

b»  = 

thickness  of  film  at  orifice 

P  = 

pressure  in  atmospheres 

X  = 

wavelength 

d  = 

radius 

SECTION  III. 

\  =  illumination 
0  =  "brightness 
N.A  =  numerical  aperture 
^  =  magnification 
d  =  depth  of  field 

Z  =  diameter  of  circle  of  confusion 
^  -  wavelength 
T  =  transmission 


appendix  A. 


As  in  Section  I. 


APPENDIX  B. 


As  in  Section  I  plus  : 

w  =  rate  of  reaction  per  unit  volume  per  second 

a  =  mass  fraction  of  reactant  A  etc. 

0=  frequency  factor  in  first  order  rate 
equation 

E.= 


energy  of  activation. 


Ref:  X.R. 465/1. 


INTRODUCTION. 


One  of  the  advantages  of  the  rocket  motor  over  other 
propulsion  units  lies  in  its  ability  to  achieve  a  high  rate 
of  conversion  of  chemical  energy  into  kinetic  energy  in  a 
relatively  small  volume.  In  order  to  maintain  the  high  rate 
and  density  of  combustion  required,  the  propellant  is  injected 
into  the  combustion  chamber  in  the  form  of  a  spray  consisting 
of  a  large  number  of  small  droplets.  This  spray  has  a  large 
surface  area  per  unit  volume  which  facilitates  the  transfer 
of  heat  from  the  combustion  gases  to  the  liquid  propellant. 

Thus  the  high  rate  of  evaporation  from  the  liquid  to  the  gaseous 
phase,  necessary  to  maintain  combustion,  is  ochived. 

In  the  case  of  multi -component  propellants  mixed  in  the 
combustion  chamber,  the  dispersion  of  each  component  is  a  step 
towards  their  mutual  interaction. 

It  is  desirable  to  have  methods  of  assessing  the  'tomising 
properties  of  injectors  for  liquid  propellants  which  can  be 
correlated  with  measurements  of  the  efficiency  of  the  injector 
in  the-  combustion  system  for  which  it  is  designed.  The  methods 
which  have  been  used  in  other  fields  for  the  measurement  of  drop 
Bize  in  sprays  are  reviewed  in  this  report,  but  first  it  is 
desirable  to  examine  the  influence  of  droplet  size  on  the 
combustion  of  propellant  sprays,  because  the  selection  of  a 
method  suitable  for  the  assessment  of  rocket  injectors  will  be 
governed  by  the  form  of  the  relation  between  drop  size  and  rate 
of  combustion.  Sprays  produced  by  methods  employed  in  com¬ 
bustion  engineering  contain  droplets  of  which  the  size  may  vary 
over  a  considerable  range  and  it  is  necessary  to  characterise 
both  the  mean  droplet  size  and  the  degree  of  deviation  from 
the  mean.  The  definition  of  these  quantities  will  depend  upon 
their  relation  to  the  mechanism  of  combustion. 

The  calculation  of  the  overall  rate  of  combustion  of  the 
propellant  Eprny  from  possible  relations  between  th-  rate  of 
combustion  of  an  individual  droplet  and  the  dimension  of  that 
droplet  is  discussed  in  the  first  section  of  this  report. 

In  order  to  do  this  a  greatly  simplified  picture  of  the  propellant 
spray  and  the  process  of  combustion  is  necessary.  It  is  assumed 
that  the  process  of  atomisation  is  instant ~ncous ,  th^t  is,  the 
propellant  is  injected  in  the  form  of  discrete  droplets,  the 
initial  size  ^nd  distribution  of  which  c  n  be  seumed,  whereas  in 
reality  the  bre  ik  up  of  the  injected  propellent  is  occurring  in 
the  combustion  chamber  simultaneously  with  its  combustion. 

It  is  to  be  expected  that  the  physical  factors  entering  into  the 
combustion  of  mono-propellants  should  be  much  simpler  than  in 
bi-propellant  systems.  Initially,  therefore,  discussion  has 
been- confined  to  the  combustion  of  mono-propellants.  The  com¬ 
bustion  of  systems  in  which  fuel  and  oxidant  are  injected  separate¬ 
ly  is  of  such  complexity  that  it  merits  a  separate  discussion 
which,  it  is  hoped,  any  form  a  subsequent  pert  of  this  dissertation. 
In  the  case  of  mono-propellants  it  is  assumed  that  each  droplet 
burns  independently  of  the  rest  of  the  spray,  and  therefore  the 
overall  rate  is  given  by  the  sum  of  the  combustion  of  all  droplets 
in  the  combustion  chamber.  It  is  then  possible  to  calculate 
the  rate  of  burning  per  drop  as  a  function  of  the  thermal  properties 
of  the  propellant  and  the  size  of  the  drop.  The  assumptions 
needed  to  give  an  algebraic  solution  ore  too  drastic  to  permit  of 

/quantitative 


quantitative  calculation  of  the  rate  of  combustion  of 
any  particular  propellant  system  but  should  give  an 
indication  of  the  dependance  of  combustion  upon  droplet 
diameter  and  therefore  permit  the  examination  of  the 
effect  of  changes  in  the  mean  droplet  size  and  inhomo¬ 
geneity  of  the  propellant  spray. 

The  process  of  atomisation  depends  not  only  on  the 
injector  but  also  on  the  conditions  existing  in  the  com¬ 
bustion  chamber  and  on  certain  properties  of  the  liquid 
propellant.  Since  it  will  be  more  convenient  to  test 
the  atomising  properties  of  injectors  in  the  absence  of 
combustion  and  with  liquids  other  than  the  propellaht  for 
which  the  injector  has  been  designed,  it  is  desirable  to 
discuss  the  effect  of  these  changes  on  the  degree  of  atom¬ 
isation,  and  therefore  the  relation  of  the  latter  to  the 
conditions  of  operation.  In  the  second  section  the 
literature  on  the  atomisation  of  liquids  is  reviewed  with 
particular  reference  to  the  effect  of  liquid  properties 
such  as  viscosity  and  surface  tension  and  to  the  operating 
conditions  such  as  injection  pressure,  air  density,  etc. 

The  third  section  reviews  methods  of  measurement  of 
drop  size  in  liquid  sprays  and  discusses  their  applicability 
to  rocket  injectors.  It  is  desirable  that  the  method  or 
methods  adopted  should  be  simple  and  speedy  and  of  general 
application  to  the  varied  types  of  spray  which  will  be 
tried  in  the  early  stages  of  rocket  development.  It  is 
unlikely,  however,  that  there  will  be  a  requirement  for  a 
purely  routine  method  of  checking  the  atomisation  of  pro¬ 
duction  types  for  control  purposes  in  the  same  way  ns  might 
be  required  for,  say,  flow  rate.  This  is  because  small" 
differences  in  the  degree  of  atomisation  will  have  n 
negligible  effect  on  the  overall  efficiency  of  the  rocket 
motor.  The  need  is  more  for  a  method  of  determining  large 
differences  resulting  from  n  change  in  type  or  scale  of  the 
injection  system  in  development  work,  and  for  research  into 
the  importance  oi  atomisation  and  injection  methods  in 
rocket  combustion.  There  is  relatively  little  information 
m  the  literature  on  this  subject,  and  whnt  there  is  lacks 
detail,  suggesting  the  need  for  further  research  in  this 
field. 
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SECTION  I. 


INFLUENCE  OF  DROPLET  SISE  ON  COMBUSTION . 


In  the  last  two  decades  a  number  of  papers  (1,  2,  3, 

4)  have  been  published  dealipg  with  the  influence  of  the 
degree  of  subdivision  of  fuel  sprays  in  compression  ignition 
engines.  It  has  generally  been  taken  for  granted  that  the 
rate  of  combustion  of  the  fuel  spray  is  directly  related  to 
the  total  surface  of  the  spray,  and  for  purposes  of  comparison 
between  different  sprays  a  mean  diameter  is  defined  such  that 
the  actual  spray  is  replaced  by  a  fictitious  homogeneous 
spray,  having  the  same  ratio  of  total  volume  to  total  surface. 
This  mean  diameter  was  first  used  by  the  German  engineer 
Souter  (l)'i  and,  to  distinguish  it  from  other  mean  diameters 
it  is  often  termed  the  Sauter  mean  diameter;  this  is  defined 
as  follows  i 


Where  ni  is  the  number  of  droplets  having  a  diameter  di, 
and  the  summations  ore  token  over  the  total  number  of  droplets 
N;  V  and  S  arc  the  total  volume  nnd  surface  respectively.  It 
is  obvious  that  an  infinite  number  of  sprays  have  the  same 
Sauter  mean  diameter,  differing  in  the  range  of  dispersion  of 
the  drop  sizes  about  this  mean,  and  attempts  were  mode  to  get 
up  some  number  characteristic  of  the  lack  of  uniformity  of  the 
spray. 

These  early  attempts  (2)  usunlly  involved  a  number  derived 
from  the  standard  deviation  of  the  droplet  distribution,  no 
methods  having  been  developed  for  measuring  the  distribution 
of  droplet  size  in  fuel  sprays  with  any  precision.  In  recent 
years,  however,  in  connection  with  the  development  of  spray 
nozzles  for  gas  turbines  (5)  and  oil  furnaces  (6),  it  has  been 
observed  experimentally  that  the  distributions  of  droplet  sizes  in 
sprays  from  swirl  type  spray  nozzles  fit  the  Roein-Rammler  type 
of  distribution  law,  first  applied  to  pulverised  coal  dust  (?) 
which  has  the  form 

R=  e'll)  <2> 


Where  R  is  the  cumulative  volume  fraction  with  diameter  • 
greater  than  d,  and  d0  and  n  are  the  size  and  distribution 
constants  respectively;  the  diameter  do  is  that  above  which 
i/e*  of  total  spray  volume  lies;  the  value  of  n  defines 

the  'spread'  of  the  distribution  and  is  found  experimentally  to 
take  values  greater  than  two,  the  larger  values  of  n  indicating 
more  uniform  distribution.  Experimentally  the  relation  has  been 
applied  to  sieve  fractionation  of  solidified  fuel  droplets  where 
relatively  large  samples  of  the  spray  can  be  examined  and  the 
volume  distribution  iB  determined  directly.  For  values  of  n 
than  two,  most  of  the  spray  volume  lies  in  the  range 
^  2,5  as  con  be  seen  from  the  table  given  on  page  35  of 

/Appendix  A 
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Appendix  A  to  this  report.  It  is  shown  there  that,  if  the 
cumulative  volume  distribution  follows  equation  (S)  then  the 
probability  distribution  functions  for  other  parameters  of  the 
spray  such  as  surface,  diameter  and  number  per  unit  volume, 
can  be  simply  expressed  in  terms  of  tabulated  functions, 
subject  to  certain  restrictions  upon  the  values  of  the  dis¬ 
tribution  constant  n. 

While  it  is  necessary  to  remember  that  an  injector  oper¬ 
ating  under  combustion  chamber  conditions  will  form  sprays  of 
markedly  different  mean  droplet  size  and  distribution  from 
those  found  in  ’cold'  tests,  it  is  plausible  to  assume  that  the 
form  of  the  distribution  law  will  be  similar.  Then  if  the 
dependence  of  the  rate  of  burning  of  the  propellant  drops  upon 
their  size  is  assumed,  the  effect  of  ch  ;ngcs  in  the  size  dis- 
constant  and  distribution  of  the  spray  can  be  examined.  For 
instance,  if  the  rate  of  burning  of  the  spray,  referred  to 
unit  volume  of  propellant,  is  proportional  to  the  surface  area 
per  unit  volume 


vC  =kS 


v-  v/ 


6k 


where  ds  is  the  S outer  mean  diameter  as  defined  inequation  (l) 
and  k  is  a  constant,  the  rate  of  burning  therefore  increases  with 
decreasing  mean  diameter.  It  can  be  shown  (P.37)  that  for  a 
spray  obeying  the  Rosin  Rommler  distribution  law 

where 


or  ns 


values  of  this  function  are  tabulated  (42) 
As  m - v  n-o  5  1  ^  i 


is  the  Gamma  function  defined  as 

foa 

-t 


rs-tt 

l(.  +  x)  =  j  x>-t 

!  (i  +x)  =  (tc'yl 


and  *n  — *  t  }  I 


As  n  varies  between  c*n  and  2,  rv«  increases  from  1  to 
approximately  2;  thus  the  rate  of  burning  of  the  spr-y  increases 
■v.'ith  decreasing  homogeneity.  This  refers  to  the  initial  rate  of 
burning  of  the  propellant  spray  and  this  as  R.P.  Probert  (a)  has 
remarked  is  of  less  importance  than  the  steady  state  of  combustion 

/in 
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in  a  chamber  into  which  propellant  is  continuously  injected. 


If  the  time  spent  by  the  propellant  in  the  combustion 
chamber  is  unrestricted  in  the  steady  state  n  constant  volume 
of  propellant  will  be  burning  with  a  volume  rote  of  burning 
equal  to  the  volume  rote  of  input.  The  droplet  size 
distribution  of  the  propellant  ond  it a  combustion  properties 
will  determine  the  volume  of  propellant  in  the  combustion 
chamber.  The  ratio  of  this  volume  to  the  combustion  chamber 
volume  represents  the  density  of  loading  of  the  combustion 
Bystem;  it  is  plausible  to  assume  that  this,  for  a  given 
propellant  and  combustion  chamber  and  injector  design,  should 
be  held  constant. 


It  is  shown  in  Appendix  A  to  this  report  that  if  the  rate 
of  burning  normal  to  the  surface  of  the  propellant  droplets 
is  independent  of  their  diameter,  so  that  the  overall  rate  of 
combustion  is  proportional  to  the  surface*  area,  then  the 
variation  of  diameter  of  any  droplet  with  time  cm  be  expressed 
as 


d  -  cl;  — kt 


(4) 


and  the  volume  of  propellant  in  the  combustion  chamber  in  the 
steady  state  is 


(5) 


assuming  that  the  state  of  the  propellant  droplet  distribution 
at  the  moment  of  injection  can  be  represented  by  (?) 

The  variation  of  this  volume  with  throughput  V  and  burning 
constant  k  is  to  be  expected.  It  is  however,  remarkable-  since 

ff7+i)  is  almost  independent  of  n  for  values  greater  than 

twom  that  it  should  depend  mainly  on  the  size  constant  d0  and 
be  but  slightly  affected  by  the  spread  of  the  distribution. 

However,  in  any  practical  case  the  time  that  can  be  spent 
by  any  droplet  in  the  combustion  chamber  is  limited  ond  drop¬ 
lets  over  a  certain  initial  size  may  be  ejected  from  the 
combustion  chamber  before  they  have  burnt  completely.  It  can 
be  shown  that  the  fraction  of  the  propellant  spray  injected  at 
time  t  =  o  remaining  unconsumLd  at  n  time  t>0  can  be 
represented  as 


(6) 


assuming  that  variation  of  drop 
represented  by 


diameter 


with 


t  ime 


can  be 


The  form  of  the  function  depends  upon  the  values  of  n  and/3 
and  cannot  in  general  be  given  explicitly.  Values  of  for 

.particular  values  of  and  n  are  best  found  from 
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equation.  (28)  (Appendix  A)  by  numerical  integration. 

Values  are  tabulated  below  for  a  range  of  values  of  n. 

It  can  be  seen  that  the  fraction  remaining  increases  with 
decreasing  n. 

Volume  fraction  of  injected  propellant  unburnt  at  time 
t  after  injection  as  a  function  of  n  nt  fcgt/ 


'll 

2 

3 

4 

Q  1 

0.10 

0.069 

0,050 

0.011 

0,(544 

0.02  19 

This  value  of  k/gl/£j3"'9  corresponds  to  complete  combustion  of 

the  homogeneous  spray  of  droplet  size  d0.  The  figures  quoted 
indicate  no  more  than  the  relatively  small  change  in  efficiency  to 
be  expected  from  changes  in  distribution  or  "spread"  of  the  drop¬ 
let  sizes  and  the  expected  conclusion: 


That  for  maximum  efficiency  with  a  limited  time  spent  in  the 
combustion  chamber  the  spray  distribution  should  be  as  homo¬ 
geneous  as  possible. 

Similar  conclusions  have  been  reached  by  R.P.  Probert  for 
the  injection  of  oil  fuels  into  gas  turbine  combustion  chambers. 
The  methods  employed  in  Appendix  A  for  the  evaluation  of  mean 
diameters  etc,,  are  generalisations  of  his  methods  for 
propellant  systems  where  the  rate  of  combustion  of  drops  can  *be 
related  to  their  diameter  by  indicial  equations  of  the  general 

f0™  _  =  k/  - 

Probert  assumes  that  the  combustion  in  gas  turbine  systems 
is  dependant  upon  the  rate  of  evaporation  of  the  fuel  which  is 
assumed  to  be  governed  by  an  equation  similar  to  the  psychro- 
metric  equation  in  its  dependence  on  linear  dimension  of  the 
drops,  i.e,  the  rate  of  evaporation  per  drop  ib  directly  pro¬ 
portional  to  the  diameter. 


Since  it  can  be  argued  that  combustion  of  propellants  in 
general  is  preceded  by  evaporation  it  is  worth  while  examining 
the  applicability  of  this  assumption  to  propellant  systems. 


The  assumption  that  the  rnte  of  combustion  of  a  drop  of 
propellant  is  proportional  to  the  first  power  of  the  drop  radius 
implies  that  the  rate  of  combustion  per  unit  surface  is  in-' 
versely  proportional  to  the  drop  radius.  Applied  to  mono- 
propellnnt  systems  this  would  seem  to  contradict  experience,  since 

n  r+ j  ’  rat o  per  unit  surface  i  s  usually 

assumed  to  be  independent  of  the  geometrical  configuration  of  the 

surface,  and  the  overall  rate  of  combustion  to  be  proportional 

resolved  Thls  ['®Enrent  contradiction  enn  be 

we  exnminG  the  process  of  combustion  in  relation  to 
the  linear  scale  of  the  combustible  drop.  It “m  bo  shorn 

/that 
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It  will  "be  shown  that  the,  two  possible  relations  for  the  rate 
of  burning  of  a  drop  of  rnonopropollnnt  /Q  *  t  or  ?. 

are  limiting  forms  of  the  relation  between  rate  of  burning  and 
linear  dimension..,.  If  one  imagines  the  liquid  drop  to  be 
surrounded  by  two  zones  then  at  the  surface  of  the  drop 
evaporation  occurs,  in  the  first  zone  the  vapour  is  heated  to 
its  ignition  temperature,  and  in  the  second,  reaction  occurs 
with  the  liberation  of  heat,  part  of  which  is  conducted  back 
to  the  surface  of  the  liquid  to  maintain  the  evaporation. 

For  convenience,  it  is  assumed  that  a  steady  state  has  been 
reached  in  which  the  heat  conducted  back  is  just  sufficient  tc 
evaporate  the  weight  of  liquid  consumed  in  the  reaction,  so  that 
the  temperature  of  the  surface  of  the  drop  is  constant. 

It  is  probable  that  such  a  steady  state  is  realised,  since  if 
the  reaction  zone  is  displaced  towards  the  liquid  surface, 
the  heat  conducted  back  will  increase,  the  velocity  'of  efflux 
of  vapour  will  increase  and  the  reaction  zone  will  move  away 
from  the  surface  towards  its  initial  position. 


It  is  assumed  that  this  constant  surface  temperature  is 
equal  to  the  boiling  point  of  the  liquid  under  the  combustion 
chamber  pressure  and  that  the  temperature  distribution  in  the 
interior  of  the  liquid  is  constant.  The  thermal  diffusivity 
of  liquids  is  of  the  order  of  io-4  cms2/ee-c  so  that  the  surface 
temperature  will  fall  rapidly  towards  the  centre  of  the  drop 
and  it  will  bo  convenient  to  include  the  hent  required  to  raise 
the  liquid  from  its  initial  temperature  to  the  surface  temper¬ 
ature  in  the  heat  of  vaporisation.  The  problem  is  then  to 
solve  the  partial  differential  e-quations  of  hent  and  moss  flow 
in  the  zone  outside  the  liquid  drop  with  certain  boundary 
conditions.  Since  we  are  interested  in  the  variation  of  rate 
of  burning  with  droplet  diameter  rather  than  in  calculating  the 
absolute  rate  of  burning,  and  since  the  form  of  the  equations 
renders  an  exact  solution  laborious,  the  problem  has  been  sim¬ 
plified  by  the  assumption  of  a  temperature  below  which  the 
reaction  is  zero.  This  temperature  is  termed  the  ignition 
temperature f  although  it  has  no  relation  to  the  temperature  of 
ignition  as  usually  measured.  Above  the  ignition  temperature 
it  .is  implicitly  assumed  that  the  variation  of  reaction  rate 
with  position  can  be  approximated  by  a  simple  function  of  the 
distance  from  the  surface,  although  this  function  is  eliminated 
irom  the  final  result  by  referring  the  rote  of  combustion  of 
the  drop  to  that  of  a  plane  surface.  If  we  than  make  the 
further  assumptions  that  the  ignition  temperature  and  the  rate 
oi  reaction  at  that  temperature  are  independentof  the  rodius  of 
curvature  of  the  propellant  surface,  and  that  parameters  such  ns 
thermal  conductivity  and  heat  capacity  can  be  assigned  a  mean 
value  independent  of  temperature,  the  following  relation  can  be 
derived  (see  Appendix  B) 


A  . 

c-m  =  •m0+- 


when  m  and  mc 
for  a  drop  of 


are  the  mass  rates  of  burning  per  unit  surface 
diameter  d  nd  a  plane  surface  respectively. 
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and  Q;= 


To  "being  the  surface  temperature  °A 

Ti  the  ignition  temperature  °A 

L  the  heat  of  evaporation  per  unit  mass 

Cp  the  mean  constant  pressure  specific 

heat  of  the  propellant  vapour  (between 
T0  and  Tj_)  per  unit  mass 

"X  is  the  mean  thermal  conductivity  of  the 
vapour  between  T0  and  Ti 

Putting  Z)'  in Q;  _  /p  it  can  be  seen  that  if 

Cp  r 

^>>^o  ^  (7a) 

and  if 

d«  -ft  *  y  +.  03  ' 1 « | 

j£ 

It  can  be  shown  that  z^0  is  the  width  of  the  zone  in 
which  the  propellant  vapour  is  heated  to  the  temperature  Tj_  in 
the  burning  of  a  plane  surface.  The  rate  of  burning  of  the 
propellant  drop  is  thus  dependent  on  the  surface  area  or  on  the 
diameter  according  to  relative  order  of  magnitude  of  the  size 
of  the  drop  and  the  width  of  the  heating  zone  for  the  case  of  a 
plane  surface. 

j£ 

The  value  of  rh0  will  depend  upon  the  properties  of  the 
particular  mono-propellant  considered;  it  will  obviously 
increase  with  decreasing  rate  of  burning.  The  latter  quantity 
can  in  theory  be  measured  experimentally  but  the  values  for 
actual  mono-propellants  are  not  known  with  any  accuracy.  For 
the  purpose  of  illustration  it  is  sufficient  to  assign  nominal 
values  to  the  rate  of  burning,  and  the  physical  constants  con¬ 
tained  in  U?  .  If  we  take  methyl  nitrate  as  nn  example  it 

is  considered  that  the  following  values  ore  plausible  at  n 
pressure  of  twenty  atmospheres. 

m0  —  1  gram/ems2  sec.  X  = 

L  =  100  cals/gram.  T0  = 

Cp  =  0.5  cals/gram0.  6,  = 


10“4  eals/ems2  sec.0  grad. 
350°A,  =  950 °A. 

3.5  and  %  &  5 u 


There  is  a  certain  amount  of  experimental  evidence  for  the  value  • 
of  fli0,  the  value  of  Ti  however  is  very  speculative.  There  is 
indirect  evidence  that  the  bulk  of  the  reaction  in  combustion 
systems  occurs  at  temperatures  close  to  the  final  temperature. 
Semenov  (9),  Wolfhnrd  and  Klnukens  (10)  Boys  and  Corner  (ll) 

£hc  t0^PGrDturo  extent  of  reaction  ns  n  function 
nnT  ?nfa?T‘StnnC&-i  T1"0!?  tht .  propellant  surface  for  the  burning  of 

^  would  appear  that  o  negligible  fraction 
of  the-  rate  determining  reaction  takes  place  below  850°A  lrnCtl0n 

Fortunately  <p  is  relatively  insensitive  to  chafes  in  the  value 

/of 
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of  Ti.  The  values  of  \  Op,  etc.,  are  probable  by  analogy 
rather  than  by  experiment ,  and  the  computed  value  of 
should  be  taken  as  an  order  of  magnitude  rather  than  a 
numerical  value.  The  inference  is  that  the  assumption  that 
the  rate  of  burning  is  directly  proportional  to  the  surface 
area  is  true  only  for  propellant  drops  above  about  200/ 
in  the  case  of  methyl  nitrate.  It  is  thus  open  to  question 
whether  this  law  can  be  applied  to  the  spray  distribution 
existing  in  rocket  combustion  chambers,  particularly  since 
although  representative  injectors  are  found  to  give  sprays 
under  atmospheric  conditions  in  which  the  bulk  of  drops  be 
between  100  and  400 ^  ,  the  distribution  under  combustion 

chamber  conditions  is  unknown.  It  is  probable  that  the 
increase  in  temperature  of  the  liquid  will  decrease  the  mean 
diameter  and  the  spread  of  the  distribution  by  decreasing  the 
liquid  surface  tension  and  viscosity,  but  the-  change  due  to 
the  alteration  in  gas  flow  conditions  cannot  be  estimated. 

If  it  is  assumed  that  the  bulk  of  the  propellant  will  have  drop 
diameters  in  a  range  between  10  and  200/*  neither  of  the 
approximate  forms  (7n)  of  the  burning  law  will  apply.  The 
form  of  equation  (7)  makes  the  calculation  of  the  burning 
characteristics  of  the  distribution  by  the  methods  outlined 
in  Appendix  B  and  referred  to  in  the  first  part  of  this  section 
somewhat  tedious.  It  can  be  shown  however,  by  a  log/log  plot 
of  m  against  d  that  the  rate  of  burning  per  drop  may  be 
expressed  with  fair  accuracy  within  a  limited  range  of  diameter 
by  a  relation  of  the  form 

M  =  kd*  (S) 

where  /?  has  a  value  between  one  and  two 

Thus  for 

I  <d  <  IO/X 
and  M  =■  1 6  ?  »o/d  <200^1 

By  the  use  of  these  approximations  the  relations  deduced  in 
Appendix  B  for  the  relevant  mean  diameter,  volume  of  propellant 
in  combustion  chamber,  volume  fraction  unconsumed  after  a 
certain  time,  etc.,  can  be  extended  to  droplet  size  ranges  in 
which  the  limiting  forms  of  the  combustion  equation  are-  inade¬ 
quate. 

So  far  we  hove  not  considered  the  possible  effect  of  any 
relative  motion  between  the  propellant  and  the  combustion 
chamber  gases.  In  order  to  achieve  n  state  of  atomisation  the 
propellant  is  given  an  initial  velocity  of  the  order  of  some 
tons  of  metres  per  second;  as  it  proceeds  along  the  chamber 
the  propellant  spray  is  decelerated  by  resistance  forces,  and 
the  combustion  chamber  gases  are  accelerated  both  by  the  shared 
momentum  of  the  spray  and  by  the  conversion  of  chemical  energy 
into  thermal  and  translational  energy.  Although  the  relation 
between  the  resistance  forces  and  the  drop  dimensions  and  relative- 
velocity  arc  known  for  conditions  under  which  no  combustion 
takes  place,  it  will  probably  be  considerably  ch  inged  for  the 
case  of  a  burning  drop,  since  the  burnt  gases  may  have  a  velocity 
normal  to  the  surface  of  similar  order  to  that  of  the  mean 
relative  velocity  between  the  drop  and  the  combustion  chamber 
gases. 
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However,  it  can  fee  concluded  that  any  relative  velocity 
will  have  "been  reduced  to  a  low  order  for  the  greater  part 
of  the  combustion  (it  will  he  noted  that  the  decelerating  force 
is  inversely  proportional  to  some  power  of  the  drop  diameter; 
therefore  this  conclusion  becomes  more  likely  the  smaller  the 
mean  drop  size  of  the  spray)  and  the  residual  relative 
velocity  will  have  only  a  small  effect  upon  the'  combustion 
of  a  mono-propellant  spray. 

This  last  conclusion  is  not,  however,  likely  for  bi- 
component  propellant  systems,  where  fuel  and  oxidant  are 
injected  separately,  and  the  rate  of  combustion  will  depend  upon 
the  rate  of  mixing  of  the  two  components.  Then  the  flow 
conditions  in  the  combustion  chamber  become  of  groat  importance. 
Since  these  arc  not  amenable  to  calculation,  a  detailed  analysis 
of  the  combustion  is  not  practicable,  and  such  conclusions  as 
can  be  drown  are  general,  and  only  approximate.  The  rate  of 
combustion  at  any  point  will  be  proportion '1  to  the  product  of 
the  concentrations  of  the  reactants  at  that  point,  and  a  stable 
zone-  of  combustion  can  be  formed  only  where  a  balanced  and 
constant  rate  of  supply  of  each  component  is  possible.  The 
latter* -will  depend  upon  the  initial  evap oration  and  subsequent 
diffusion  and  convection  of  each  component  from  the  liquid 
drops  to  the  combustion  zone,  and  it  is  pertinent  to  examine 
how  these  processes  depend  upon  the  diameters  of  the  individual 
droplets. 

Since  the  equations  determining  heat  flow  and  mass  flow 
are  similar  in  form,  the  requirement  for  a  simultaneous  solution 
of  these  equations  will  not  affect  the  order  of  the  dependence 
upon  the  dimensions  of  the  system,  and  only  one  need  be  considered. 
The  rate  of  mass  flow  across  a  spherical  surface  *oncentric 
with  the  drop,  within  which  the  rate  of  reaction  is  inappreciable, 
will  be  directly  proportional  to  the  radius  of  the  drop  only  if 
the-  distance  of  the  surface  from  the  drop  is  large  compared  to 
the  drop  radius;  thus  the  linear  dependence  of  evaporation 
upon  drop  size  assumed  by  Probert  will  hold  only  if  reaction 
occurs  at  a  distance  from  the  drop  largo  compared  to  the  drop 
radius;  ns  the  reaction  zone  distance  approaches  the  drop 
diameter,  the  rote-  of  evaporation  will  become  proportional  to  the 
surface  area  of  the  drop. 

Under  conditions  in  which  the  drop  has  n  relative  trans¬ 
lational  motion  with  respect  to  the  gaseous  medium,  the  rate  of 
evaporation  is  increased  by  n  term  which  is  dependent  upon  the 
Reynolds  number  of  the  flow  conditions. 

Prom  measurements  of  the  rate  of  evaporation  from  spheres 
in  a  flowing  medium  by  Prozling  (is) ,  Kramers  (13)  has  obtained 
the  following  relationship  between  the  dimensionless  numbers 
representing  the  ratejof  evaporation,  the  flow  conditions  and 
the  properties  of  the  gaseous  mediums 

^  =  2  +  055  (ft)5  (Re)4  (9) 


where  m  is  the  rate  of  evaporation  per  unit  surface  ner  unit 
concentration  difference  aSd  d  the  Srop  diameter? D  She  Effusion 

/coefficient 
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coefficient.  The  Prandl  number  Pr,  representing  the 
properties  of  the  gaseous  medium,  is  defined  as  ;  RrVgp 

*  A 


\being  the  thermal  conductivity  and  the  kinematic 
viscosity,  and  the  Reynolds  number  representing  the  condition 
of  flow  as  „  , 

Re  = 


being  the  relative  flow  velocity. 


The  constant  term  in  equation  (9)  is  equivalent  to  the 
psychrometric  equation  where  the  overall  rate  of  evaporation  is 
proportional  to  the  drop  diameter;  as  the  Reynolds  number 
becomes  large  the  second  term  will  contribute  a  term  varying 
with  the  three  halves  power  of  the  drop  diameter.  Since  for 
most  gases  and  vapours  the  Prandl  number  is  of  the  order  of 
unity,  the  second  term  will  be  appreciable  for  Reynolds  numbers 
greater  than  one.  For  kinematic  viscosities  of  the  order  of 
2  x  10”1  e.g.s.  and  drop  diameters  of  200  u  this  condition  will 
arise  if  the  relative  velocity  is  greater 'Than  one  metre  per 
second.  Thus  it  is  unlikely  that  under  combustion  chamber 
conditions  a  linear  dependence  of  rote  of  burning  upon  droplet  * 
diameter  will  hold. 


It  is  doubtful  whether  it  is  possible  for  multicomponent 
propellant  systems  in  general  to  go  further  than  the  statement 
that  the  mte  of  burning  will  increase,  ind  the  time  for  com¬ 
bustion  decrease  with  decreasing  mean  drop  size  of  the  components. 
Apparent  exceptions  may  be  found  even  to  this  generalisation, 
for  instance  with  reacting  bicomponent  systems  injection  in  the 
*  form  of  a  spray  could  be  disadvantageous  compared  with  simple 
impingement  of  the  two  liquids  if  the  energy  released  by  partial 
reaction  has  a  more  powerful  dispersing  force  than  any  physical 
means. 


/SECTION  II. 


15 


SECTION  II. 


THTT;  AT  Or  ~Tn  AT'ION  OF  LIQUIDS:  A  REVIE':/ . 


The  influence  of  the  physical  properties  of  fuels  upon 
the  state  of  atomisation  produced  by  various  injection  methods 
is  important,  both  from  the  standpoint  of  injector  design,  and 
for  the  testing  and  comparison  of  different  individual  systems. 

In  many  coses  it  will  not  be  desiroble  in  "cold"  tests  to 
use  the  propellant  for  which  the  injector  wos  designed,  but  to 
use  n  liquid  which  is  similar  to  the  propellant  in  the  relevant 
physical  properties  or,  to  adjust  other  more  eoeily  variable 
parameters  in  such  a  way  ns  to  give  n  dynamically  similar  system. 

In  the  past  twenty  years  a  considerable  amount  of  research 
has  been  done  on  the  effect  of  design,  operating  conditions  and 
the  physical  constants  of  the  sprayed  liquid  on  spray  character¬ 
istics,  notably  by  Lee  (14)  on  injectors  for  Diesel  engines, 

Hottel  ( IS )  on  spray  nozzles  for  oil  furnaces,  Nukiyano  and 
Tnmasavn  (16)  on  nir  atomisers,  and  in  England  by  Shell,  Lucas 
and  N.G.T.E.  on  injectors  for  gas  turbines.  A  comprehensive 
review  of  fuel  spray  phenomena  relating  to  internal  combustion 
engines  has  been  made  by  Muraezcw  (17)  and  an  extensive  list  of 
references  is  given  by  Roesch  and  Rose  (13). 

However,  in  spite  of  all  this  it  is  clear  that  the  precise 
relations  between  the  degree  of  atomisation  and  the  possible 
variables  are  still  unknown;  in  view  of  the  apparent  complexity 
of  the  process,  and  the  technical  difficulties  of  isolating  one 
variable  from  another,  thiE  is  not  surprising.  The  necessity 
for  the  statistical  treatment  of  results  is  apparent,  since  the 
size  of  any  particular  drop  is  not  uniquely  determined,  but 
governed  by  a  probability  relation;  this  must  be  borne  in  mind 
when  considering  the  relationships  that  have  been  put  forward 
as  governing  ntomisation,  and  explains  the  fact  that,  while  there 
is  general  agreement  on  the  direction  of  the  variation  of  mean 
droplet  size  with  the  various  parameters,  there  is  very  little 
on  the  magnitude  of  the  variation.  Further,  many  of  the 
result s  quoted  in  the  literature  ore  applicable  only  for  the 
conditions  and  techniques  employed  in  the  original  research. 

It  is,  nevertheless,  agreed  that  no  fundamental  difference 
exists  in  the  mechanism  of  the  three  most  common  methods 
employed  for  the  atomisation  of  liquids,  namely,  'solid1  or 
high  pressure  injection,  centrifugal  injection,  and  nir  or  gns 
atomisation.  These  differ  chiefly  in  their  field  of  applic¬ 
ability,  and  the  method  used  for  initial  dispersion  of  the 
liquid.  The  forces  causing  ntomisation  ore  the1  inertia  of 
the  liquid  and  the  resistance  forces  between  the  liquid  and 
gaseous  medium.  The  first  process  relies  on  turbulence  to 
spread  the  liquid  into  high  velocity  sheets,  so  that  the  resist¬ 
ance  forces  may  complete  the  process  of  ntomisation,  and  requires 
high  injection  pressures  to  give  the  high  flow  velocities  re¬ 
quired;  the  second,  as  its  nomenclature  would  imply,  uses 
rotational  forces  to  achieve  the  initial  dispersion  of  the  jet, 
and  therefore  operates  with  only  moderate  pressure  heads;  the 
last  method  imparts  the  high  velocity  to  the  gaseous  medium,  and 
needs  only  a  small  liquid  pressure  head;  it  is  most  useful  for 
low  rates  of  liquid  input,  and  for  combustion  systems  containing 
one  gaseous  component.  With  correct  design  and  choice  of 
operating  conditions,  nil  three  methods  will  trive  equal  degrees  of 
ntomisation. 
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/The 


The  theory  of  the  break  up  of  liquid  jets  has  been 
examined  by  Rayleigh  ( 19)  who  concluded  that  for  a  column  of 
liquid  in  a  state  of  unstable  equilibrium  due  to  the  action 
of  capillary  forces,  the  wave  length  of  the  disturbance  lead¬ 
ing  most  rapidly  to  disintegration  is  given  by  =  4.508  x  the 
diameter.  This  was  confirmed  by  observations  on  low  velocity 
liquid  jets  which  showed  a  rotationally  symmetric  periodic 
variation  in  diameter,  increasing  with  increasing  distance 
from  the  orifice,  until  individual  drops  were  formed  by  rupture 
at  the  constructions;  the  wave  length  of  the  disturbance,  and 
therefore  the  diameter  of  the  drops  formed  wse  found  to  be 
proportional  to  the  jet  radius. 


By  the  use  of  spark  photography,  Haenlein  (3)  discovered 
that  at  higher  jet  velocities  the  periodic  contraction  of  the 
jet  became  rotationally  asymmetric,  and  at  still  higher 
velocities,  was  replaced  by  a  wave-like  jet  of  constant  cross- 
sect  ions,  which  tended  to  break  at  the  crests  and  the  troughs 
through  the  action  of  air  resistance  in  increasing  the  amplitude 
of  the  wave;  with  increase  of  velocity,  n  critical  value  was 
reached,  .at  which  the  disruption  of  the  jet  into  drops  appeared 
to  occur  at  the  orifice.  Castleman  ( 31 )  explained  the  latter 
process  in  terms  of  the  Rayleigh  theory  of  the  instability  of 
liquid  jets  by  assuming  that  small  portions  of  the  jet  surface 
were  drawn  out  by  air  resistance  forces  into  -fine  ligaments, 
which  then  broke  off  from  the  main  surface  and  contracted  by 
surface  tension  forces  into  small  drops.  This  was  contrary  to 
current  theory,  which  supposed  that  high  pressure  injectors 
relied  on  the  internal  turbulence,  and  the  suddenly  released 
pressure  energy  of  the  liquid,  to  effect  atomisation,  so  that 
the  mechanism  was  essentially  different  from  that  of  nir  atom¬ 
isation,  and  would  occur  even  in  vacuo. 


The  experiments  of  Lee  (22)  on  spark  photomicrography  of 
diesel  engine  sprays  showed  that  Casteman's  views  on  atomisation 
were  essentially  correct;  he  observed  that,  nt  low  (one 
thousandth  of  atmospheric)  nir  densities,  the  jet  is  spread  out 
into  sheets  by  orifice  turbulence,  the  influence  of  orifice 
conditions  and  liquid  viscosity  ond  velocity  becoming  very  marked* 
but  that  the  sheets  collapse  under  the  notion  of  surface  tension  ’ 
forces  to  give  n  relatively  coarse  spray. 


As  the  air  density  is  increased,  the  stages  in  drop  form¬ 
ation,  proposed  by  Castleman,  are  observed;  the  distance  from 
the  orifice  within  which  a  stable  state  of  atomisation  is 
reached  decreases;  at  normal  air  densities  the  effect  of 
increasing  density  is  more  on  this  distance  than  on  the  final 
degree  of  atomisation,  and  at  high  densities  the  effect  of 
rifice  design  and  fuel  viscosity  is  small,  except  in  so  far  as 
they  may  doorcase  the  initial  jet  velocity. 


OT1,  Studied  the  spray  formation,  both  by  direct  photography, 

.nd  by  the  collection,  measurement  and  counting  of  the  droolet^ 

£rn?  +h the-°r^fiCV  ^eign,  the  liquid'' in je ct ion  pressure 
head,  the  chamber  air  density  and  the  viscosity  and  surface 

tension  of  the  liquid  spray.  His  conclusions  were  summarised 

"  ,hCn  VOri?^cs  nrc  held  constant,  the  degree 

isintegr.it ion  of  the  jet  (n)  increases  with  the  distance  ' 
from  the  nozzle,  until  the  disintegrating  forces  due  to  the 

suffiMrnrtnCi  y  tCtwe+£  the  nir  nnd  thc  fuel  are  no  longer 
sufficient  to  overcome  the  resisting  forces  due  to  the  surface 

tension  and  viscosity  of  the  fuel,  fb)  increases  with  increase 
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of  air  density;  (c)  increases  with  increase  of  jet  velocity; 
(d)  decreases  with  increase  of  orifice  diameter;  (e) 
decreases  with  increase  of  fuel  viscosity;  (f)  decreases  with 
increase  of  fuel  surface  tension;  (g)  increases  with  increase 
of  fuel  turbulence.  Fuel  turbulence  accelerates  the  dis¬ 
integration  of  the  jet  by  ruffling  its  surface  close  tcbthe 
orifice,  but  has  relatively  little  disintegrating  effect  in 
itself" . 


Lee,  probably  wisely,  did  not  attempt  to  formulate  any 
precise  relationship  between  a  mean  drop size  and  the  other 
parameters;  such  attempts  as  have  been  made  by  other  workers 
are  based  upon  a  somewhat  crude  dimensional  analysis. 

It  is  assumed  that  there  are  two  forces,  the  inertia  of  the 
liquid  and  air  resistance,  tending  to  break  up  the  jet,  while 
viscous  and  surface  forces  resist  disintegration.  These 
forces  may  be  expressed  dimensionally  in  terms  of  in  velocity, 
d  representative  length,  the  surface  tension,  viscosity  and 
density  of  the  liquid  and  density  of  the  gaseous  medium 
(the  effect  of  the  viscosity  of  the  gaseous  medium  is  assumed 
to  be  negligible).  Then; 


(1) 

(2) 

(3) 

(4) 


A  '  j  Q 

the  inertia  force  given  by  fi  xr  has 
dimensions  & 

the  air  resistance  force  given  by  kpA-? 
dimensions  £ 

the  viscosity  force  given  by  ,u  A  jt-  has 
dimensions  v  f  "L. 


the  surface  tension  force  given  by  yl- 
dimensions  y| 


ha  s 


has 


In  the  orifice,  and  at  low  air  densities,  terms  (l)  and 
(3)  may  be  taken  to  be  predominant,  and  the  condition  for 
dynamic  similarity  is  that  their  ratio  is  constant 


Condition  A. 


L/Ig 

V 


In  the  orifice  L  con  be  identified  ns  the  orifice  radius  and 
condition  A  implies  that  the  Reynolds  number  of  the  flow 
should  be  held  constant,  i.e.  an  increase  in  fuel  viscosity 
can  be  compensated  by  increasing  the  flow  velocity,  maintaining 
the  same  typo  of  flow  and  degree  of  turbulence  in  the  orifice."’ 
In  the  free  jet,  L  can  be  taken  to  be  the  drop  diameter  and 
it  has  been  assumed  that,  for  viscous  liquids  at  low  air 
densities,  the  most  probable  drop  diameter  is  directly  pro¬ 
portional  to  the  viscosity,  and  inversely  proportions  to  the 
jet  velocity. 


At  higher  air  densities,  terms  (2)  and  (4)  v/ill  predominate, 
and  therefore  we  obtain 


Condition  B.  _ £_ 

This  implies  that  the  most  probable  drop 
proportional  to  the  surface  tension,  and 


size  will  be  directly 
inversely  proportional 
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to  the  air  density  and  the  square  of  the  jet  velocity. 

It  should  ho  remembered  that  ratio  B  is  constant  only  when 
the  Reynolds  number  with  respect  to  the  gaseous  medium  is 
high,  and  the  resistance  of  the  medium  is  independent  of  the 
Reynolds  number.  Treibnigg  (23)  considered  this  condition 
to  hold  for  the  maximum  possible  drop  size  since,  for 
stability,  the  average  air  pressure  over  the  surface  of  the 
drop  must  be  balanced  by  the  surface  tension  pressure. 

Scheubel  (24)  modified  this  condition  by  introducing  the  ratio 
between  surface  tension  and  viscous  forces,  i.e,  dU  <p(X  \ 

but  did  not  check  the  relationship  between  drop  size  r' *  ' 

and  viscosity  experimentally. 


Sass  (85)  found  that  the  reciprocal  relationship  between 
drop  si2G  and  air  density  was  approximately  true  for  high  air 
densities,  but  that  there  was  a  smaller  rate  of  increase  of 
drop  size  with  decreasing  density  at  atmospheric  densities. 

Lee  (26)  and  Re tel  (27)  found  that  the  mean  drop  size  goes 
through  n  minimum  value  as  the  air  density  is  increased,  but 
this  conclusion  has  been  criticised  on  the  grounds  that,  in 
Retel's  experiment,  the  increase  in  drop  size  may  be  due  to 
decreasing  jet  velocity,  since  he  was  operating  with  a  fixed 
injection  pressure,  and  in  Lee's  case  the  high  spray  density 
obtained  may  have  favoured  re-combination  of  the  drops. 

Sauter  and  Scheubel  found  the  drop  size  in  air  atomisers  to  be 
approximately  proportional  to  the  surface  tension,  and  inversely 
proportional  to  the  square  of  the  air  velocity.  Nukiynmo  and 

Tamasava  (loc.cit)  deduced  the  following  empirical  formula  i 

1 

as = f )  +  Mvjjy  .  («ga.)  (»> 

where  Qq  and  are  the  volume  flow  rates  of  liquid  and  air 
respectively. 


Lewis  and  Richards  (28)  have  concluded^  that  Sauter* s  date 
do  not  conflict  with  this  equation.  They  also  find  that, 
although  as  it  stands  the  empirical  equation  of  Nukiyama  and 
Tamasava  is  not  dimensionally  homogeneous,  with  suitable 
assumptions  to  obtain  the  appropriate  value  of  Qa ,  it  gives 
the  correct  order  of  magnitude  of  drop  size  when  applied  to 
pressure  and  centrifugal  injectors. 


Merrington  and  Richardson  (29)  studying  the  break  up  of 
betn  stationary  jetg  and  jets  discharged  from  nozzles  having  a 
velocity  with  respect  to  the  air,  found  that  at  high  air/liquid 
velocities  the  volume  medium  drop  size  at  a  large  distance1 
f rom  tin;  orifice  woe  independent  of  the  orifice'  diameter  and 

thetmeB!nftheT,™2S1$  ?£  wos  proportional  to 

when  break  up  occurred  during  the  time  of  fall,  S  P  ’ 
With  pressure  and  centrifugal  injectors,  the  liquid/nir 


-  19 


/velpcity 


velocity  is  not  measured  directly,  but  must  be  calculated 
from  the  liquid  density  and  the  pressure  drop  across  the 
nozzle  by  means  of  Bernoulli's  equation 

a.AF*  (n) 


£ 


In  the  case  of  viscous  liquids  the  actual  velocity  will  be 
what  less  than  that  given  by  this  equation,  owing  to  the 
dissipation  of  energy  in  the  orifice.  For  plain  orifices 
is  possible  to  calculate  the  velocity  from  the  volume  flow 


some- 

it 

rate 


Q_=  k.Area,  jaAjF^ 

where  k  is  the  discharge  coefficient,  and  is  tho  product  of  the 
velocity  coefficient,  i.c.  the  ratio  of  the  actual  velocity  to 
that  given  by  Bernoulli's  equation,  and  the  contraction  co¬ 
efficient,  which  is  the  ratio  of  the  jet  area  to  the  orifice 
area.  At  Reynold's  numbers  greater  than  the  critical  value 
(S,OOC  -  3,000)  the  contraction  coefficient  can  be  taken  to  be 
close  to  unity  and  the  velocity  coefficient  equal  to  the 
measured  discharge  coefficient.  The  latter  has  a  usually 
constant  value  of  between  0.7  -  0.9  above  the  critical  Reynolds 
number,  the  actual  value  depending  upon  the  orifice  length/ 
diameter  ratio  and  fluid  viscosity. 


With  centrifugal  nozzles  the  contraction  coefficient  is 
usually  less  than  unity,  due  to  the  presence  of  an  air  core  in 
the  orifice,  and  is  a  function  of  the  viscosity.  It  has  been 
shown  by  Watson  (30)  that,  neglecting  viscosity,  the  discharge 
coefficient  can  be  calculated  from  the  ratio  of  tangental  to 
axial  velocity  imparted  to  the  liquid;  and  therefore  from  the 
design  of  the  nozzle.  Thus  it  might  be  possible  to  estimate  the 
velocity  losses  from  the  ratio  of  the  actual  to  the  colculotod 
discharge  coefficient . 

However,  Taylor  (31)  has  shown  that,  even  for  liquids  of 
kinematic  viscosities  of  the  order  of  that  of  water,  the  major 
part  of  the  liquid  flow  forms  n  boundary  layer,  and  it  would 
appear  from  his  calculations  that  constant  discharge  co-efficients 
and  similarity  of  flow  will  only  obtain  when  the  Reynolds  number 
is  kept  constant.  In  view  of  this,  it  is  not  easy  to  compare 
the  results  of  investigators  on  centrifugal  nozzles,  where  tho 
drop  size  is  represented  ns  a  function  of  the  pressure  drop 
across  the  nozzle  with  the  equations  relating  drop  size  to 
relative  velocity  and  other  parameters.  Hottel  and  Longwell,* 
using  centrifugal  nozzles  having  cone  angles  between  60°  and  120° 
with  burning  oils  having  viscosities  between  0.1  and  0.9  stokes 
at  injection  pressures  between  50  and  300  pounds  per  square  inch, 
found  that  their  results  could  be  expressed  as: 

c^>stn%  e6,'rosV 

~ =  ap™*  as) 


do  being  the  mean  drop  size  defined  by  the  Rosin  Rnmmclin 
equation 


R= 


r0  the  radius  of  orifice  and  oc  is  the  spray  cone  angle.  The 
latter  is  itself  dependent  upon  the  viscosity,  decreasing  by 
about  20‘S  over  the  range  of  viscosities  studied.  The  introduction 
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of  the  term  Sin  %  which  helps  to  linearise  the  dependence 
of  drop  diameter  on  orifice  radius  is  justified  "by  Longwell 
on  the  grounds  that  the  thickness  of  the  conical  sheet  of 
liquid,  formed  at  a  distance  1  from  the  orifice  along  its 
surface,  is  proportional  to  'TfetoO  “  J$ 


where  t0  is  the  thickness  at  the  orifice.  Then,  if  it  is 
assumed  that  the  size  of  the  drops  formed  by  "break  up  of  the 
sheet  is  proportional  to  its  thickness,  and  if  to  is  very  much 
less  than  ro,  i.e.  if  the  discharge  coefficient  is  small,  then 


i.e.  inversely  as  the  ratio  of  the  length  of  the  liquid  sheet 
to  its  thickness  at  the  orifice.  This  ratio  is  assumed  "by 
Longwell  to  "be  independent  of  the  scale  and  design  of  the 
nozzle,  and  to  be  a  function  only  of  the  infection  pressure  and 
liquid  viscosity.  However,  if  the  initial  thickness  is 
independent  of  pressure,  as  it  should  be  if  the  discharge  co¬ 
efficient  is  constant,  the  1  :ngth  of  the  conical  sheet  should 
increase  with  increasing  pressure  to  effect  a  decrease  in  drop 
size.  This  is  not  borne  out  by  observation,  for  whereas  the 
pressure  increases,  the  point  of  break  up  of  the  sheet 
approaches  the  orifice.  Work  in  England  has  led  to  the  con¬ 
clusion  that  the  drop  size  is  a  minimum  when  the  cone  angle  is 
about  90°  and  is  not  linearly  proportional  to  the  orifice  size. 
Needham  (33)  haB  concluded  that,  for  a  number  of  injectors  of 
this  optimum  angle,  the  drop  size  can  be  expressed  by 


Q'25. 
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(13) 


where  the  pressure  range  covered  was  6  to  125  pounds  per  square 
inch ,  _ 

B  /2.  A.F  ' 

since  Q,_=  kirVM-£- 

(14) 
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This  apparent  divergence  of  opinion  on  the  dependence  of 
drop  size  on  orifice  size  and  injection  pressure  is  no  doubt 
due  in  part  to  different  methods  of  treating  the  data,  but  it 
is  indicative  of  the  difficulty  of  formulating  a  generalised 
theory  of  atomisers  of  even  one  type,  the  centrifugal  mozzle. 

For  the  testing  and  comparison  of  rocket  injectors,  it  is 
clear  that  there  is  not  yet  enough  information  available  on 
the  effect  of  the  physical  properties  of  liquids  on  drop  size 
to  permit  the  simulation  of  one  liquid  by  another,  allowing  for 
the  change  in  one  property  by  adjustment  of  another.  If,  for 
reasons  of  convenience,  it  is  desired  to  use  a  substitute  liquid 
each  relevant  physical  property  should  be  reproduced  os  closely 
as  possible.  In  so  far  ns  the  properties  of  the  gaseous  medium 
into  which  the  liquid  is  sprayed  are  concerned  simulation  is  not 
possible.  In  the  combustion  chamber  the  gas  density  may  be 
from  two  to  four  times  that  of  air  at  atmospheric  pressure. 
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The  mean  gas  temperature  may  lie  between  1000  and  2000 °A, 
therefore  the  gas  viscosity  will  be  some  four  to  eight  times 
larger.  The  changes  in  the  kinematic  viscosity  will  be 
rather  less  owing  to  the  increased  density.  The  effect  of 
the  increased  temperature  of  the  liquid  on  its  viscosity  and 
more  particularly  on  its  surface  tension  must  be  considerable. 
The  surface  temperature  of  the  liquid  will  be  equal  or  near 
to  the  boiling  point  under  the  combustion  chamber  pressure; 
this  may  approach  the  critical  temperature  at  which  the  surface 
tension  vanishes. 

The  flow  condition  in  the  gaseous  medium  in  a  combustion 
chamber  will  differ  considerably  from  those  of  an  unignited 
spray  in  the  open  air.  It  has  been  observed  that  the  shape 
and  distribution  of  sprays  are  modified  by  the  shape  and  size 
of  the  chamber.  Further,  under  conditions  of  forced  evap¬ 
oration  the  velocity  of  flow  normal  to  the  surface  of  the 
liquid  mny  be  comparable  to  the  tangential  velocity  between  the 
liquid  surface  and  the  gaseous  medium;  this  will  reduce  the 
resistance  forces  tending  to  break  up  the  drop. 

It  is  not  possible  to  assess  the  overall  effect  of  these 
conditions  on  the  atomising  properties  of  nn  injector  and 
there  seems  to  be  n  wide  field  for  experimental  research  into 
the  functioning  of  injectors  under  conditions  approaching  those 
in  the  combustion  chamber. 
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SUCTION  III 


THE  MEASUREMENT  OF  SPRAY  DROPLET  SIZE:  A  REVIEW . 


It  Is  obvious  from  the  preceding  discussion  that,  while 
knowledge  of  the  drop  size  distribution  of  injector  sprays 
is  desirable  in  so  far  as  size  distribution  must  affect  the 
efficiency  and  extent  of  the  combustion,  and  while  it  will 
be  profitable  to  attempt  to  correlate  this  information  with 
proof stand  tests,  a  high  degree  of  accuracy  is  not  required. 

On  the  other  hand  in  view  of  the  uncertainty  as  to  the  order 
of  relationship  (quantitatively)  between  drop  size  and  com¬ 
bustion  characteristics  it  will  be  desirable  to  measure  the 
size  distribution  rather  than  any  mean  parameter  such  as  specific 
surface. 

Bearing  these  facts  in  mind,  it  is  useful  to  review  briefly 
the  methods,  both  direct  and  indirect,  which  have  been  used  to 
measure  drop  sizes  in  sprays.  In  the  former  the  droplets  in 
a  representative  sample  of  the  spray  are  caught,  measured  and 
counted,  while  indirect  methods  rely  on  the  measurement  of  a 
property  of  the  spray  which  can  be  related  to  some  mean  dimension 
The  advantage  of  direct  methods  is  the  ease  with  which 
distribution  curves  can  be  obtained  from  the  data.  The  dis¬ 
advantages  are  the  excessive  time  and  labour  required  for  the 
assessment  of  the  spray,  technical  difficulties  in  catching  the 
spray  in  a  form  suitable  for  measurement,  and  the  sampling 
problems  which  arise  for  non-homogene ous  sprays. 

On  the  other  hand  indirect  methods  in  general  do  not  give 
ony  information  on  the  dispersion  of  drop  sizes  about  the  mean 
value,  require  considerable  ingenuity  for  application  and  -ore 
liable  to  somewhat  large  instrumental  errors;  but  results 
can  be  obtained  without  the  expenditure  of  excessive  time  and 
labour  after  the  initial  setting  up  of  the  apparatus. 

*  The  direct  methods  of  spray  droplet  size  determination  used 
by  previous  investigator's  (33)  differ  mainly  in  the  means 
adopted  to  catch  samples  of  the  spray  in  a  form  suitable  for 
subsequent  counting  and  measurement  of  the  individual  droplets. 
The  simplest  method  is  to  expose  a  glass  microscope  slide  to 
the  spray  for  such  a  time  as  to  give  a  convenient  density. 

If  an  uncoated  slide  is  used,  the  liquid  drops  wet  the  glass 
surface  forming  piano  convex  lenses  of  which  the  volume  can  be 
estimated  from  their  diameter  and  focal  length.  This  method 
is  tedious  and  can  only  be  successfully  applied  to  non-volatile 
liquids;  further  the  droplets  must  not  exceed  a  certain  size 
or  they  will  tend  to  break  up  on  impact  with  the  slide.  A 
simple  method  which  can  be  applied  to  large  drops  is  to  catch 
them  on  sheets  of  absorbent  paper,  the  liquid  is  preferably  dyed 
so  that  the  stains,  which  can  be  previously  calibrated  using 
drops  of  known  size,  are  easily  observed  (29). 

Slides  coated  with  a  layer  of  viscous  liquid  in  which  the 
droplets  embed  themselves,  retaining  their  spherical  shape,  hev> 
been  used  for  both  oily  and*  aque-ous  sprays,  gelatinised 
glycerine  and  mixtures  of  paraffin  jelly  and  paraffin  oil 
respectively  being  employed.  If  the  medium  is  previously 
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saturated  with  the  sprayed  liquid  and  the  slide  covered  with 
a  cover  slip  after  exposure,  the  records  ore  stable  over  a 
period  of  several  days  and  even  the  smallest  droplets  can  he 
observed. 

Surfaces  coated  with  carbon  and  magnesium  oxide  smoke  will 
show  identations  due  to  the  impact  of  droplets  which  may  be 
related  to  the  size  of  the  droplets  by  previous  calibration. 
Slides  coated  with  a  soluble  dye  have  been  shown  to  be  sensitive 
to  droplets  of  the  order  of  1  /j  ,  Dobie  (34)  has  used  dishes 
containing  castor  oil  diluted  with  alcohol  to  give  a  density 
similar  to  that  of  the  sprayed  liquid  and  has  also  caught  the 
droplets  on  photographic  plates  coated  with  an  oil  and  vaseline 
mixture  obtaining  a  permanent  record  by  subsequent  exposure  to 
light,  removal  of  the  oily  film  with  ether  and  development  of 
the  light  sensitive  emulsion. 

These  methods  need  some  form  of  mechanical  shutter  between 
the  spray  and  the  slide  if  a  density,  small  enough  to  prevent 
coalescence  of  the  drops  upon  the  slide,  is  to  be  obtained. 

With  dense  sprays  and  short  times  of  exposure  there  will  be  n 
danger  of  break  up  of  drops  upon  the  edges  of  the  shutter,  and 
with  droplet  diameters  of  the  ranges  existing  in  injector  sprays 
there  is  a  high  probability  of  break  up  of  the  larger  drops  on 
impact  with  the  elide.  Although  the  method  lends  itself  to 
the  simultaneous  determination  of  variations  of  spray  density 
and  droplet  distribution  at  various  planes  in  the’ spray  volume 
it  is  not  considered  that  it  is  of  general  application  for 
injector  testing. 


An  ingenious  method  of  catching  the  spray  in  a  form  suit¬ 
able  for  either  direct  assessment  or  for  the  indirect 
determination  of  surface  area  has  been  used  by  Longwell  in 
America  and  by  the  Shell  Laboratories  in  England,  '  The  former 
has  collected  oil  sprays  in  n  freezing  bath  "of  solid  carbon 
dioxide  dissolved  in  alcohol;  the  oil  droplets  solidify  and 
the  particle  size  distribution  is  obtained  by  fractional  sieving; 
the  use  of  low  temperatures  has  obvious  disadvantages  on  the 
experimental  side.  Joyce  (5)  has  examined  the  sprays  produced 
by  gas  turbine  injectors  by  substituting  for  the  fuel  used 
operationally,  a  suitable  gr^’de  of  paraffin  wax  heated  to  such 
a  temperature  that  it  has  physical  properties  similar  to  the 
fuel.  The  liquid  wax  spray  solidifies  in  flight  and  the  wax 
particles  are  collected  in  water.  They  may  then  be  fractionally 
sieved,  or  measured  and  counted  on  a  microscope  slide  (35). 
Successful  attempts  have  been  made  to  suspend  the  particles  in  a 
suitable  liquid  and  determine  the  specific  surface  by  measuring 
the  optical  transmission  of  the  suspension  and  it  is  possible 
that  the  size  distribution  curve  could  be  obtained  bv  sedi¬ 
mentation  (36). 


This  method  has  the  considerable  advantage  that  the  sample 
obtained  is  representative  of  the  whole  spray,  and  with  the 
sieving  method  can  be  made  much  larger  than  is  possible  with 
methods  involving  counting.  The  validity  of  the  assumption 
that  the  solid  paraffin  wax  spray  is  identical  with  that  given 
by  liquid  oil  is  open  to  question:  it  will  depend  on  the  rate 
of  cooling  ot  the  spray  and  the  time  taken  for  a  stable  state  of 
atomisation  to  be  reached.  Hollow  wax  droplets  have  been 
observed  in  wax  sprays,  due  probably  to  the  slight  solubility 
of  the  pressurising  gas  in  the  liquid  wax,  and  it  would  seoni 
desirable  that  a  check  against  a  method  of  known  accuracy,  using 
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"both  liquid  wax  and  the  operational  liquid,  should  be  made 
before  the  method  can  be  accepted  as  reliable.  It  is  not 
very  suitable  for  sprays  of  high  droplet  density  where  partial 
coalescence  of  the  wax  spheres  tend  to  occur,  and  it  is 
difficult  to  see  how  this  method  can  be  applied  to  injectors 
designed  to  operate  with  liquids  resembling  water  in  their 
physical  properties,  or  for  impinging  jet  atomisers. 

The  difficulties  inherent  in  the  application  of  the  above 
methods  to  injector ■ Bprays  have  led  to  the  study  of  the 
possibilities  of  photography  of  sprays  in  flight,  using  high 
intensity  short  duration  illumination  to  “freeze"  the  droplets 
in  space.  Although  spark  photography  has  long  been  used  to 
obtain  information  on  the  mechanism  of  atomisation,  it  does 
not  seem  to  have  been  much  used  for  the  measurement  of  droplet 
sizes,  probably  because  the  demands  for  high  resolution  and 
large  depth  of  field  are  irreconciliable.  However,  the 
advantages  in  flexibility  of  application  of  this  method  are 
considerable.  In  distinction  to  most  other  methods,  there  is 
no  interaction  between  the  subject  and  the  method  of  observation’ 
it  can  be  applied  to  any  point  in  the  spray;  in  rocket  injectors 
considerable  interest  may  be  centred  on  the  unstable  states  of 
dispersion  close  to  the  orifice  or  at  the  point  of  impingement 
of  jets.  Providing  transparent  windows  can  be  used," it  can  be 
applied  to  atomisation  under  conditions  other  than  atmospheric, 
and  in  high. velocity  gas  streams.  Simultaneously  with  drop 
size  determination  it  will  give  information  on  the  spray 
density  at  any  point  and  with  slight  modification  on  the  droplet 
velocities.  The  main  non -technic U  disadvantage  of  the  method 
ns  applied  to  injector  assessment  lies  in  the  tendency  to  focus 
attention  on  the  individual  droplet  at  n  certain  point  in  space 
^nd  time,  whereas  the  time  average  over  the  total  spray  volume 
is  more  relevent.  This  necessitates  care  in  sampling  inhomo¬ 
geneous  or  intermittent  sprays. 


The  experimental  difficulties  lie  in  the  use  of  ultra  short 
exposure  times  to  prevent  blurring  of  the  image  of  the  drop  by 
its  motion  relative  to  the  camera  and  obtaining  adequate  reso¬ 
lution  without  excessive  limitation  of  the  depth  of  field,  The 
first  is. not  a  serious  problem  at  the  present  time;  for  a  drop 
moving  with  a  velocity  of  t on  metre s  per  s econd  an  exposure  of 
1  micro  second  will  limit  the  movement  during  exposure  to  10 
micron;  since  a  limit  of  resolution  of  the  order  of  20  micron 
is  adequate,  the  exposure  time  of  one  or  two  microseconds.  Riven 
by  standard  'pArditron'  equipment,  should  be  satisfactory  fo? 
most  sprays*  J 

.  The  sprays  encountered  in  injector  research  contain  droplets 
having  diameters  of  between  50  and  1000  micron;  the  weight 
fraction  of  spray  outside  this  range  is  small  and  for  any 
particular  spray  the  bulk  is  contained  in  a  ten  or  twenty  fold 
dianieter  range.  Thus  for  convenience  in  measurement  a  ten  to 
litty  fold  linear  magnification  will  be  required. 

?cvfral  f masons  it  is  desirable  to  perform  this  magni- 
thC  in  tK°  st^s>  using  a  magnification  of  about  two  on 

the  photographic  plate  and  obtaining  the  desired  degree  of 

SageS!f;elaLdU?r^r^Pr0JtCJ10^  Th£  U^?nf?ion°of  the 

image  is  related  to  brightness  of  the  source  B,  the  numerical 
*Pby  Se  ?elaUoSi:,e0tlVe  K'A;  ana  the  Sesree  of  ™S»ifioation 
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T'nc  level  of  nllurainntion  required  is  determined  by  the  speed 
of  the  emulsion  used,  rind  to  obtain  high  contrast  this  is  kept 
low.  In  view  of  the  short  exposures  needed  the  source _ bright¬ 
ness  used  is  the  maximum  obtainable.  In  order  to  obtain 
maximum  depth  of  field  the  numerical  aperture  is  made  small 
subject  to  the  limitation  of  diffr  ction  effects;  thus  it  is 
necessary  to  restrict  the  initial  degree  of  magnification.  A 
lower  limit  for  this  is  fixed  by  the  resolution  -'nd  grain  size 
of  the  emulsion,  the  former  for  process  emulsions  is  about 
ten  M  ,  and  although  the  required  degree  of  resolution  is  unlike¬ 
ly  to  bo  less  than  twenty//  it  is  undesirable  to  have  an  initial 
magnification  of  less  than  unity.  Further,  the  projection 
distance  becomes  inconveniently  long  if  second  stage  magnifica¬ 
tions  of  more  than  fifty  are  required. 

The  depth  of  focus  is  determined  by  the  degree  of  resolu¬ 
tion  required;  the  indeterminacy  in  the  measurement  of  the 
diameter  of  a  spray  droplet  will  be  approximately  equal  to  the 
resolving  power,  i.e.  the  least  distance  between  two  object 
points  nt  which  they  may  be  distinguished.  For  a  spray  of 
weight  median  diameter  of  two  hundred  micron,  a  resolution  of 
twenty  micron  is  sufficient.  The  depth  of  field  is  related 
to  the  circle  of  confusion  z,  and  the  numerical  aperture  N.A. 
by  the  relation: 

d  (16) 

while  a  relation  exists  between  zl ,  the  limit  of  resolution, 

N.A.  and  the  wave  length  of  the  illumination  X  ,  representing 
diffraction  by  the  aperture: 

Z  *  i  (&A) 


There  is  no  object  in  using 
than  Msl,  therefore,  eliminating 

d  = 

■A 

Thus,  if  zl  is  put  equal  to  20  ju 
depth  of  focus  will  be  1.6  m.m, 
lution  is  obtained  only  over  a  relatively  small  depth  of  the 
spray.  This  depth  is  independent  of  the  optical  parameters  of 
the  lens  system,  and  is  dependent  only  on  the  wave  length  of  the 
illumination.  It  is  impracticable  to  decrease  this  other 
than  by  using  a  blue  sensitive  emulsion  having  n  maximum  sensi¬ 
tivity  nt  0.45  jjl  . 

In  view  of  the  lack  of  information  on  the  practicability  of 
the  photographic  method,  experimental  work  has  been  started  at 
E.R.D.E.  During  early  experiments  an  attempt  was  made  to  limit 
the  depth  of  field  photographed  by  restricted  illumination  at 
right  angles  to  the  camera  axis.  The  droplets  were  recorded  by 
the  light  reflected  into  the  camera  lens.  However,  no  short 
duration  source  of  sufficient  brightness  was  found,  and  it  was 
decided  to  concentrate  work  on  axial  illuminating  systems,  the 
drops  being  silhouetted  against  the  light  field.  The  use  of 
slits ^to  limit  the  depth  of  spray  in  the  camera  field  was 
examined,  but  although  it  was  found  possible  to  prevent  break-up 

/of 
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a  circle  of  confusion  z  lt-ss 
N.A.  from  equations  16  and  17, 

(18) 

and  X  -  0.5//  the  maximum 
Thus  a  given 'degree  of  reeo- 
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of  the  spray  on  the  edges  of  the  slits  by  making  them  of 
absorbent  material  and  limiting  the  time  of  exposure  to  the 
spray  by  means  of  a  shutter,  it  was  decided  that  the  introduction 
of  mechanical  devices  into  the  spray  would  remove  the  special 
advantages  of  the  photographic  method. 


When  a  droplet  field  of  a  depth  large  compared  to  the  depth 
of  field  of  the  camera  system  is  photographed,  two  sources  of 
error  are  introduced.  Outside  the  focus  zone  the  definition 
falls  off  gradually,  so  that  in  the  assessment  of  the  spray 
photograph  the  operator  has  to  judge  the  point  at  which  measure¬ 
ment  becomes  impracticable.  As  a  result,  the  large  droplet 
sizes  ore  estimated  over  n  greater  depth  of  the  spray  than  the 
smaller  droplets  and  the  small  droplet  fractions  are  under- 
weighted.  For  a  given  error  in  estimation  of  droplet  diameter, 
the  fractional  volume  error  will  be  proportional  to  the  recip¬ 
rocal  of  the  droplet  diameter.  On  the  other  hand,  there  is 
a  tendency  to  under-estimate  the  diameter  of  drops  out  of  focus 
the  error  increasing  with  lack  of  sharpness  of  the  image. 


The  distortion  introduced  into  the  distribution  curve,  due 
to  these  effects,  was  examined  by  photographing  a  suspension  of 
paraffin  wax  spheres  in  water  contained  in  parallel  sided  glass 
absorption  colls  of  1  find  2  cm,  thickness.  The  distribution 
curveof  the  wax  suspension  had  been  previously  determined  by 
counting  and  measuring  a  sample  taken  on  n  microscope  slide. 

1  itli  the  cell  of  1  cm,  thickness,  the  distribution  curve  obtained 
was  not  significantly  different  from  that  given  by  the  slide 
count.  With  the  2  cm.  cell  the  apparent  weight  median  diameter 
had  decreased  by  approximately  20 M  >  and  the  slope  of  the 
distribution  curve  showed  a  slight  increase.  It  was  concluded 
that  the  accuracy  of  the  method  was  sufficient  for  sprays  of 
size  constant  greater  than  200 yU  .  With  sprays  of  higher 
dispersion  the  method  is  comparative  rather  than  absolute. 

*he  assessment  of  the  photographic  record  of  the  spray  is  best 
performed  by  projecting  the  image  on  to  a  ground  glass  screen 
The  counting  and  measurement  of  a  sufficient  number  of  drop  * 
images  is  somewhat  laborious  since  more  than  five  hundred  must 
be  counted  in  order  to  obtain  reproducible  results.  It  would 
+  hAP°SSi+^e  demise  a  photoelectric  scanning  device  operating  on 

the  negative  which  would  sort  and  count  the  drops  in  a  number  of 
size. ranges  with  sufficient  accuracy,  but  the  time  and  labour 
required  to  design  and  construct  this  would  be  prohibitively 
large.  Manual  measurement  and  counting  can  be  speeded  up  bv 
the  use  of  an  adjustable  comparison  scale,  the  movement  of  which 
a  mechanically  connected  to  a  multi-position  switch  selecting 
the  appropriate  counting  register.  A  counting  aid  operating  on 
this  principle  is  described  in  Part  II  of  this  report;  g  on 

thr  r^o?l0u:Ld+be  P°s?ible»  by  measuring  the  mean  transmission  of 
the  negative,  to  obtain  a  value  for  the  total  surface  area  of 

the  spray  sample  photographed.  The  problems  arising  will  be 

f  the  review  0f  0th£r  raothods  lnair.0?  Bseessment 


Indirect  methods  rely  on  the  measurement  of  some  property  of 
the  spray  which  can  be  related  to  a  mean  diameter  Th  LthL 
which  would  seem  to  be  most  applicable  ?o  inject ^  sprays  is  th, 
measurement  of  the  optical  transmission  of  the  spray.  Then 

jjr  ®  fpray  ^ns^ty  and  path  lenSth  such  that  over-lopping  of  the 
droplets  in  the  direction  of  the  optical  axle  can  he£og??=tcd- 

T=  |-i  v£=  54 

23s  (19) 
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as Burning  the  droplets  are  opaque.  Here  S/V  is  the  surface 
area  per  unit  volume  of  liquid,  v  is  the  volume  concentration 
of  liquid,  and  1  is  the  optical  path  length.  If  T  is  less 
than  0.8,  the  error  due  to  overlapping  of  the  droplets  cannot 
he  neglected,  and  the  exponential  equation 


(20) 


must  he  used.  This  holds  only  over  path  lengths  within 
which  the  values  of  v  and  ds  are  constant.  Thus  it  con  he 
applied  only  to  homogeneously  distributed  sprays;  otherwise 
the  overall  transmission  must  he  greater  than  0.8.  This 
limits  the  possible  applicability  of  the  method  to  injector 
sprays  to  those  of  low  density  or  spatial  homogeneity. 

With  values  of  the  transmission  coefficient  greater  than  0.9 
the  error  in  ds  for  a  given  error  in  the  transmission 
coefficient  becomes  considerable  (it  is  nine  times  the  error 
in  transmission  coefficient  at  T  =  0.9).  The  measurement  of 
the  ratio  of  the  volume  of  liquid  to  the  volume  of  spray 
presents  some  difficulties.  Souter  (37) ,  using  the  method 
to  investigate  air  atomisers,  assumed  the  volume  concentration 
of  liquid  to  be  equal  to  the  ratio  of  the  volume  flow  rates 
of  liquid  and  air,  thereby  assuming  the  velocity  of  the  air 
and  the  liquid  to  be  identical. 

In  the  case  of  pressure  and  centrifugal  injectors,  the 
quantity  most  susceptible  to  measurement  is  the  volume  flow 
rate  across  unit  surface  at  any  point  in  the  spray.  To 
derive  the  concentration  from  the  flux  density,  the  mean 
velocity  of  the  spray  normal  to  the  surface  must  be  known  cr 
assumed.  This  necessity  for  measurement  of  spray  velocities 
renders  the  method  inconvenient  for  the  routine  testing  of 
pressure  and  swirl  injectors. 


Sauter  also  used  the  charge  carried  by  drops  emitted 
from  o  nozzle  held  at  a  fixed  electrostatic:  potential  to 
measure  uhe  moan  diameter  with  respect  to  surface;  the  experi™ 
mental  difficulties  make  this  method  of  historical  interest 
only.  A  method  which  might  be  employed  as  a  rough  check 
upon  the  diameter  and  homogeneity  of  the  sprnye  relies  on  the 
diffraction  of  light  from  a  point  source  observed  through  the 
spray.  If  monochromatic  light  is  used,  and  the  spray  is 
homogeneous,  n  series  of  concentric  diffraction  rings 
surrounding  the  source  can  be  observed.  The  half  angle 
subtended  at  the  eye  by  the  nth  dark  ring  is  related  {38)  to 
the  drop  diameter  d  and  the  wave  length  ;v  by  the  relation: 

ci  -  _ 

‘  ’  Shv&"  (21) 

The  definition  of  the  rings  deteriorates  with  increasing 
inhomogeneity. 
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SUMMARY. 


The  text  is  divided  into  three  sections,  the  first  of 
which  is  n  discussion  of  the  relntion  of  propellant  atom¬ 
isation  to  combustion  in  rocket  motors,  the  second  is  a  review 
of  the  literature  on  the  atomisation  of  liquids  and  the  third  . 
n  review  of  methods  of  measuring  drop  size  in  sprays.  The 
reason  for  this  order  of  presentation  is  that  the  method  even¬ 
tually  selected  for  assessing  the  dispersive  power  of  liquid 
propellant  injectors  will  he  determined  by  the  importance  of 
atomisation  in  the  combustion  process  of  the  rocket  motor  and 
by,,  the  variation  of  droplet  size  with  changes  in  the  physical 
properties  of  the  propellant  and  the  external  conditions. 

A  spray  consisting  of  a  number  of  spherical  droplets  of 
different  sizes  can  be  described  by  two  quantities,  one 
representing  a  mean  droplet  diameter  and  the  other  the  'spread* 
of  droplet  size  obout  this  mean  diameter.  It  has  been  found 
experimentally  that  the  distribution  of  droplet  size  in  the 
sprays  given  by  some  types  of  atomisers  can  be  represented  by  the 
Rosin-Ronmler  expression; 

R=  e.i) 

This  equation  defines  the  'size  constant',  do,  and  the  dis¬ 
tribution  constant,  n.  It  is  shown  that  if  a  relation  between 
the  rate  of  combustion  of  n  droplet  of  propell-nt  and  its 
diameter  can  be  found  the  combustion  characteristics  of  a  spray 
can  be  defined  by  do,  "h  ,  and  the  droplet  diameter  constants 
in  the  rate  of  combustion  equation.  By  analogy  with  the 
combustion*  of  solid  propellants  it  might  be  thought  that  the 
rate  of  combustion  of  a  droplet  of  a  monoprop ell ant  would  be 
proportional  to  the  surface  area,  that  is  to  the  square  of  the 
diameter  of  the  droplet.  This  conclusion  is  shown  to  be  true 
only  when  the  radius  of  curvature  of  the  droplet  surface  is 
small  compared  with  the  dimensions  of  the  preheating  zone 
which  must  exist  in  the  vapour  phase  surrounding  the  droplet. 

In  theory  the  variation  of  rate  of  combustion  can  be  calculated 
from  the  physical-chemical  constants  of  the  propellant  system; 
in  practice  the  calculation  is  laborious  and  the  value  of  the 
constants  of  any  particular  system  uncertain.  An  approximate 
treatment  xe  attempted  which  indicates  that  for  droplets  of 
the  order  of  diameter  to  be  expected  in  the  combustion  of 
liquid  propellants  in  rocket  systems  the  rate  of  combustion 
per  drop  will  be  proportional  to  some  power  of  the  diameter 
between  one  and  two. 

i.e.  M=  ?  |  </3<2 

The  value  of  can  be  estimated  if  the  rate  of  burning  per  unit 
area  of  plane  surface  of  the  monopropellnnt  is  known.  As 
this  increases  A  tends  to  its  limiting  value  of  two. 

In  the  case  of  bipropellant  systems  it  is  impossible  to 
make  any  estimate  of  the  value  of/3  without  n  detailed  examination 
of  the  mechanism  of  the  process  of  combustion  for  particular 
systems.  The  influence  of  changes  in  the  state  of  gas  flow  in  s 
the  combustion  chamber  will  make  the  experimental  study  of  the 
influence  of  droplet  size  more  difficult  than  in  monopropellnnt 

/  systems. 
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syst ems. 


Assuming  that  the  rate  of  combustion  per  drop  is  given 
by  p'-i  =  kjkcr  and  that  the  distribution  of  droplet  sizes  in 

the  spray  can  be  represented  by  the  initial 


rate  of  combustion  of  the  spray  per  unit  volume  of  propellant 


is  given  by 


_,AT„ 


M  -  Kg  i  17,1  3 z§\ 

v  ^  hi  / 


3 -0 


This  expression  indicates  that  the  initial  rate  of  combustion 
of  the  spray  increases  with  decreasing  size  constant  and 
decreasing  homogeneity  (decreasing  values  of  n). 


Under  steady  operating  conditions  the  overall  rate  of 
combustion  must  be  proportional  to  the  rate  of  input  of 
propellant  (and  equal  to  the  latter  for  maximum  efficiency). 
The  volume  of  propellant  in  the  combustion  chamber  is  given  by 


K  a  being  defined  by  the  equation  for  the  variation 
diameter  with  time: 


* 


3-/3 

5=  d;  -  lf/8t 


of  droplet 


Vc c  increases  with  increasing  d0  indicating  that  if  the  density 
of  combustion  is  kept  constant  the  volume  of  the  combustion 
chamber  must  increase  with  increasing  dG;  on  the  other  hand  it 
is  not  very  dependant  on  the  value  of  n.  The  fraction  of 

propellant  unburnt  at  a  time  t  after  injection  is  .mor,e.o in¬ 
dependent  upon  the  inhomogeneity  of  the  spray.  It  can  be 
expressed  as  a  function. of  the  parameter  ket/d^  At  a 

value  of  t  which  makes  this  expression  unity'  a  change  of  n  from 
4  to  2  increases  the  fraction  of  propellant  unburnt ~ by  a  factor 
of  2  for  «  |  and  by  a  factor  of  5'^  for  -  2.  . 

This  formal  treatment  indicates  the  direction  of  the  effect 
of  changes  in  the  fineness  and  inhomogeneity  of  the  propellant 
spray.  The  uncertainty  in  the  appropriate  value  of b the  index 
in  the  rate  of  burning  equation  and  the  simplifying  assumptions 
used  in  this  analysis  make  more  precise  deductions  fruitless. 


The  next  section  on  the  process  of  atomising  reveals  a 
similar  uncertainty  in  our  knowledge.  It  is  clearly  time  that 
the  fineness  of  a  spray  increases  with  decreasing  density,  surface 
tension,  and  viscosity  of  the  propellant,  and  decreases  with 
increase  of  orifice  size  and  decrease  of  injection  pressure  drop 
or  combustion  chamber  gas  density.  It  is  difficult  to  assess  the 
exact  state  of  the  atmosphere  in  which  the  atomisation  takes 
place  in  a  rocket  motor  and  the  effect  of  this  state  upon  the 
process  of  atomisation.  The  temperature  and  flow  conditions  are 
quite-  unlike-  those  encountered  in  tests  in  which  the  spray  is  not 
ignited.  It  is  probable  that  the  dispersion  is  assisted  by  the 
temp g ro ture  and  that  the  resultant  spray  is  finer  than  that 
given  by  the  _ same  injector  in  ’cold*  tests.  The  process  of 

18  not  ins‘5ontnrieoue»  ondf  since  it  takes  place  during 
ustion,  may  never  be  complete.  It  seems  desirable  to 

/investigate 
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investigate  this  process  under  combustion  conditions. 

Clearly  a  high  degree  of  precision  will  not  he  needed  in 
the  assessment  of  the  atomising  properties  of  rocket  injectors. 

The  requirement  is  for  a  method,  flexible  in  application,  which 
will  give  at  least  a  semi-quant i to tive  measure  of  the  distribution 
of  droplet  sizes  and  which  can  be  applied  to  the  incipient  stages 
of  atomisation  as  well  as  to  the  fully  developed  spray.  For 
these  reasons  after  reviewing  the  methods  which  have  been  used 
for  droplet  size  measurement  we  have  concentrated  on  photographic 
means  of  recording  the  spray  in  flight.  It  is  shown  that  this 
method  is  capable  of  greater  accuracy  than  is  octunlly  required 
both  in  the  estimation  of  the  diameter  of  an  individual  droplet 
and  in  the  measurement  of  the  distribution  curve  of  the  spray  as 
n  whole.  It  is  particularly  suited  to  the  Btudy  of  the  unstable 
state  of  the  spray  close  to  the  injector  orifice  and  of  atom¬ 
isation  under  conditions  of  high  temperature,  high  air  density 
etc.  The  theory  of  the  photography  of  sprays  is  discussed  in 
detail.  The  technique  and  apparatus  developed  at  E.R.D.E.  is 
described  in  Port  II  with  the  experimental  results  obtained  for 
representative  injectors.  It  is  intended  to  correlate  teets  on 
injectors  giving  widely  varying  mean  droplet  diameters  with 
proofstand  measurements  on  monopropellants. 

The  majority  of  the  problems  raised  have  not  been  solved. 

It  is  hoped  however,  that  it  will  serve  to  stimulate  interest  and 
further  research. 

The  author  wishes  to  thank  Mr.  L.A.  Wiseman  for  many  helpful 
discussions  and  Miss  M.  Hooper  for  assistance  with  the  numerical 
calculations. 
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APPENDIX  A 


THE  USE  OP  THE  ROSIN-RAMBLER  DISTRIBUTION  LAW 
IN  THE  REPRESENTATION  OF  PROPELLANT 

SPRAYS. 

In  dealing  with,  size  distribution  of  uniformly  shaped 
particles- one  is  usually  interested  in  the  volume  fraction 
of  the  distribution  having  dimensions  between  any  two 
arbitrary  limits  rather  than  in  the  number  of  particles  in 
this  size  range;  particularly  so  when  the  form  of  the  dis¬ 
tribution  is  obtained  by  sieving  methods  and  the  method  of 
presentation  of  distribution  curves  is  a  plot  of  the  fraction 
cr  percentage  of  the  volume  having  a  diameter  greater  than  a 
length  d,  against  values  of  d.  From  observations  on 
powdered  coal  Rosin  and  Rammler  noted  that  the  distribution 
curves  could  be  made  linear  by  plotting  the  double  logarithm  of 
the  reciprocal  of  the  volume  fraction  oversize,  R,  against  the 
logarithm  of  the  diameter  d.  This  indicates  that  the 
distribution  curve  can  be  expressed  as  : 


or,  log, log  5  =  l0S  b  +  n  log  d  +  log. log  e 

n  being  the  slope  of  the  straight  line  obtained.  Since  n 
defines  the  ’spread'  of  the  points  it  is  termed  the  distribution 
constant.  It  is  possible  to  define  the  constant  b  in  terms 
of  the  volume  median  diameter,  the  diameter  above  which  lies 
one  half  the  total  volume,  but  it-  is ' mathematically  neater  to 
define  a  diameter  do  by  the  equation: 

S  =  e"[£)  n  (1) 

Then  l/etl1  of  the  volume  has  a  diameter  greater  than  do. 

This  expression  has  been  found  to  fit  the  droplet  distributions 
existing  in  centrifugal  sprays,  and  as  cited  (P,10)  has  been 
used  in  conjunction  with  a  relation  between  the  rote  of  evapor¬ 
ation  and  the  diameter  of  volatile  droplets  to  deduce  certain 
conclusions  about  the  burning  of  oil  sprays.  The  distribution 
constant,  n,  which  defines  the  slope  of  the  distribution  curve 
has  been  found  to  take  values  between  two  and  four  for  oil 
sprays  produced  by  swirl  atomisers.  It  is  useful  to  examine 
the  effect  of  variation  of  n  on  the  volume  and  numerical  dis¬ 
tribution.  The  variations  of  the  volume  friction  oversize  at 
various  arbitrary  values  of  d/do  are  tabulated  below: 


d/do 


n 

= 

2 

0.1 

0.990 

0.  2 
0.961 

0.4 

0.852 

0.8 

0.527 

1.2 

0.237 

1.6 

0.077 

2.  C 
0.018 

2.4 

0.003 

n 

= 

2.5 

0.997 

0.982 

0.905 

0.570 

0.206 

0.C39 

0.0C4 

0.000 

n 

- 

3.0 

0.999 

0.992 

C .  938 

0.599 

0.177 

0,007 

0.000 

0.000 

n 

= 

3.5 

1.00U 

•  .997 

0.961 

0.640 

0.144 

C  .006 

0.000 

0.000 

n 

4.C 

1.CC0 

0.998 

0.975 

0.664 

0.126 

0.001 

C.CO'"' 

o.coo 

Thus  for  sprays  having  distributions  of  this  type,  the  volume 

/fraction 
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fraction  having  diameters  outside  the  ranges  0.1- d0/  d  2.4d0 
for  n  =  2  and  0.2do  <(  d  <("  1.6d0  for  n  =  4  is  negligible. 

Putting  d/d0  =  X  and  differentiating  equation  (l)  with 
respect  to  X,  we  obtain  _  . 

dR  =  -nX"  e-X  dX  (8) 

dR  is  the  volume  fraction  between  X  and  X  +  dX.  Putting  this 
equal  to  dX ,  'V*  being  the  fractional  volume  distribution 
function  and  differentiating  again 

^  ---nX^'V*  (n-t-'iOT)  (3) 

Equating  (3)  to  zero,  it  follows  that  the  maximum  volume  of  the 
spray  occurs  at  a  diameter  given  by  ; 

X=t#  <*> 

It  can  be  seen  that  the  volume  distribution  curve  has  a 
maximum  at  a  real  value  of  X  provided  that  n  >  1.  Prom  equation 
(2)  the  same  condition  holds  forVxto  vanish  at  X  =  0;  this 
implies  that  the  distribution  law  (l)  applies  to  real  distributions 
only  for  values  of  n  greater  than  one. 


The  distribution  functions  of  other  parameters  of  the  system 
can  be  derived  from  equation  (2)  e.g. 


Sv  the  surface  distribution  function  = 


'vVx 

«  X 


$vaax 


Lx  the  diameter  distribution  function 


6V* 


Nx  the  number  distribution  function 


6  Vx_ 

IT  NX1  X 


■J 


Sv,  Lv,  Nv  being  the  total  surface,  diameter  and  number  respect¬ 
ively  per  unit  volume.  The  integral  in  the  denominator  of  each 
expression  can  be  evaluated  by  its  relation  to  the  gamma  function 


which  is  identical  with  the  generalised  factorial 
function  (x) I  Values  of  this  function  can  be  found  in  Jahnke- 
Emde  (42) . 

Y-n-‘  -x 

Since  from  equation  (2)  =■  _Ti  ^  ® 


then 


\  dX  = 


&0  ^ 
-nX  e“  dX 
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It  can  be  shown  that 


fvVxdX  =  (lr1!  >  x>" 


it  follows  that 


7x"'“*V*7X  =  (V-')! 


J 


provided  that  n^oC 


By  differentiation  of  /X  with  respect  to  X  it  can  he  shown  that 
it  has  a  maximum  at  a  value  of  X  given  hy  X  =  | Ttr-j-* - 1 yfr 


When  the  distribution  function  corresponding  to  the  value 

of  «  has  no  maximum  for  real  values  of  X  and  is  finite  at  X  =  0. 
When Ti ^ G<  the  distribution  function  is  infinite  at  X  =  0  and 
the  definite  integral  (8)  cannot  be  evaluated.  If  the  total 
value  of  the  parameter  i.e.  S,  L,  N  etc.  in  (5,6,7)  or  cumulative 
values  are  required  for  *n<c*+i  it  is  necessary  to  impose  a 
lower  limit  on  X.  Then  in  the  range  <x  <ti  ^ot+i  the  integral 
can  be  evaluated  from  tables  of  the  incomplete  gamma  function 
and  for  n'  =  oC  from  values  of  the  exponential  integral;  for 
I'  numerical  integration  is  necessary.  If  the  rate  of 

combustion  etc.,  per  drop  is  given  by  ketr  ,  the  rate  of  com¬ 
bustion  per  unit  volume  is  '  then  in  the 

notation  of  P.3,  total  rate  over  all  droplet  sizes  is 


M  =  K, 


Tt 


■  M 
\ 
/ 


V, 

6^ 


(11) 


and  the  mean  diameter  for  this  dependence  of  rate  upon 

diameter  is  defined  by  the  equation  : 


Overall  rate  of  combustion 


=  t>  kj<f  vx-'a  j 


thus 


(11a) 


(12) 


I 


When  ft  =  1,  this  represents  the  mean  diameter  considered  by 
Probert  and  when  fi  =  2,  it  is  identical  with  the  Snuter  mean 
diameter. 


Prom  equation  (9) 


7  =  id< 


7-  V) 


h 


(13) 


when  the  summation  is  replaced  by  integration  over  all  possible 
values  of  droplet  diameter  ranging  from  zero  to  infinity.  As 
"h  — ^  3-/3  the  denominator  of  equation  13  and  therefore  the 

rate  of  combustion  per  unit  volume  tends  to  infinity,  and  it  is 
necessary  to  assume  a  lower  limit  of  droplet  size. 
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If  ft  =  1 


-.-4 


(14) 


we  can  compare  the  function  with  that  derived  from 


(15) 


for  values  of  n  "between  2  and  3,  and  X(  an  arbitrary  limit* 


If  we  put  y  =  Xn  in  the  integral  of  equation  15,  it 
"becomes 


(16) 


or  if  v<yz  j 


(17) 


This  expression,  the  incomplete  factorial  function,  can  "be 
found  from  tables;  these  give  values  of 

then  equation  17  is  best  formulated  as 


when  n  =  2,  the  integral  reduces  to  the  exponential  integral 


values  of  which  are  tabulated* 

In  the  calculation  of  cumulative  volume  curves  from  the 
measurements  of  the  number  of  drops  having  diameters  between 
arbitrary  limits,  the  error  due  to  the  finite  width  of  the 
diameter  increment  may  be  appreciable  when  the  diameter  range 
is  divided  into  less  than  ten  increments.  The  cumulative 
fraction  oversize  is  calculated  from  the-equntion: 


(18) 


Nj_  being  the  number  of  drops  having  diameters  between  di  and 
dq  “  A  dq.  The  mean  diameter  ^jis  taken  to  be  the  arithmetic 
moan  of  the  limiting  values  when  t^d^  is  app  reel  able  compared  to 
di*  The  error  due  to  this  assumption  has  been  examined  by 
comparing  the  computed  curves  with  an  original  cumulative  curve 
obeying  the  Rosin-Rnmmler  formula. 
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Prom  equation  (7)  the  fraction  of  the  total  number  of 
droplets  between  diameters di  and  dg 

IS  [  udx/fv* dX  X1  bol,lg  &  oto>  as) 

Jx  ~WJ  If 

A)  O 

Prom  equations  (2)  and  (9)  this  is  equal  to 

(l-^,X2)l -(l'^>Xl)!y(l-3'nV-  i  11  >3 

This  can  be  written  as  ICm^-IKp)  * 

0  ~=~  and  p  =  -5. 

Given  the  values  of  Xq  etc.,  and  n,  this  can  be  evaluated 
from  tables  of  the  incomplete  garoiin  function..  * 


Taking  n  =  4, 

so  that 

p  =  -0*75  and 

u  =  zX* 

* 

X 

* 

0.1 

0.3 

0.5 

0.7 

0.9 

1.1 

1.3 

1.5 

1.7 

u 

0.002 

0.016(2) 

0.125 

0.480 

1.31(2) 

2. 92(8) 

5.712 

10.12(5) 

16.704 

Ky.fi 

0.11034  0.33044 

0.54485 

0.73747 

0.89143 

0.96658 

0.99404 

0.99954 

0.99999 

M(u,P) 

0.11034  0.22010 

0.21441 

0.19261 

0.15396 

0.08172 

0.02913 

0.00550 

0.00045 

ax3 

- 

0.008 

0.064 

0.216 

O.512 

1.0 

1.728 

2.744 

4.096 

Av 

- 

0.00618 

0.04816 

0.14601 

0.  276o6 

0. 28682 

O.17069 

0.05300 

0.00649 

R 

y4 

1.0 

0.9933 

0.9456 

0.7997 

0.5230 

0.2362 

0.0595 

O.OO65 

0.0000 

ex 

0.9999 

0.9919 

0.9394 

0.73ou 

0.5189 

0.2322 

0.0578 

0.0061 

0.0002 

*b*\0A 

1 

19 

162 

131 

61 

40 

17 

4 

2 

The  last  line  gives  the  difference  between  values  of  the  cumulative  volume 
fraction  derived  from 


and  the  values  computed  for  finite  intervals  of  width  0. 2X  using  the 
arithmetic  mean  of  the  interval  for  the  calculation  of  the  volume  fraction. 

If  we  consider  the  injection  of  a  propellant  spray  into  a  combustion 
chamber 
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*  See  ’'Tables  of  the  Incomplete  Oamma  Function" .  Karl  Pearson, 
Stationery  Office,  1922, 


of  unlimited  size  at  a  rate  V  volumes  of  propellant  per  second, 
a  steady  state  will  be  reached  in  which  the  volume  and  droplet 
size  distribution  will  be  such  that  the  rate  of  combustion  is 
balanced  by  the  rate  of  input.  This  volume  and  distribution 
will  depend  upon  the  burning  characteristics  of  the  propellant 
and  can  be  calculated  from  this  and  a  knowledge  of  the 
distribution  of  the  spray  at  the  moment  of  injection. 


If  the  mass  rate  of  combustion  per  drop  is  proportional  to 
a  power  /$  of  the  diameter  then  the  rate  of  change  of  diameter 
with  time  is  preport ional  to  d  ,  By  integration  it  follows 

that  if  di  is  the  diameter  at  t  =  o  the  diameter  at  time  t  is 
given  by  g 

d'  8  =  d;3  i  l^t 

where  ^  is  a  constant  ( 20 ) 

and  provided  that  ^  ^-3 

Then  the  ratio  of  the  volume  of  spray  made  up  of  drops 
initial  diameter  d’f  at  time  t  to  its  initial  volume  is 


(21) 


.  If  M  is  volume  distribution  function  of  the  injected  spray 
at  d  =  di  then  V- VpS^.St  is  the  volume  between  di  and  di  -  6d 
injected  in  time  St  and  the  volume  in  the  combustion  chamber  due 
to  the  injection  of  drops  size  of  this  range  is 


o 


The  integration  is  performed  from  the  time  of  injection  to  t,=  ^ 
the  time  at  which  the  drops  vanish.  J 


If  the  distribution  of  the  injected  spray  follows  the  Rosin- 
Rnmmler  expression 


R  = 


-x 


then  the  volume  fraction  distribution 

v'0*'  -X11 
is  -tiA.  G 


where  X  - 


and  the  volume  of  propellant  in  the  combustion  chamber  from 
droplet  sizes  between  X  and  X  +  dX 

m  yJ:8  y3'e  ■v:  dx 

W.,  A 


and  the  total  volume  from  the  whole  range  of  spray  sizes 

ii  vd»/- 

“  114-  i _ 

T\ 


S  -/3 


(23) 
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/The 


-he  volume  fraction  distribution  function  of  the  propellant 
in  the  chamber  is 


(24) 


In  the  two  limiting  cases  of  propellant  combustion  &  =  1 
and  2  then  the  volume  of  propellant  in  the  combustion  chamber 


will  be  given  by  _ 

V  =  |  ^  V) 

(35) 

on4 

(26) 

Neither  of  these  gamma  functions  is  very  sensitive  to  changes  of 
n  above  the  value  two, 

A  quantity  which  is  important  when  the  time  which  the 
propellant  spray  has  for  combustion  is  limited  is  the  volume 
fraction  of  the  propellant  injected  at  t  =  o  remaining  at  time  t. 


From  equation  21  this  is  given  by 


putting  pf* 

zL  S.V  vSsfc)  „  -wr 

<J.  A  -ar 


this  equation  for\rt  is  a  function  of 


It  con  be  evaluated  by  numerical  integration. 


(2V) 


(83) 


and  of  n  and jQ 


Footnote. 


This  formula  differs  in  the  value  of  the  numerical  constant 
from  that  given  by  Probert  (loc.cit).  His  derivation  differs 
in  the  diameter  variable  chosen  and  in  the  paths  of  integration; 
the  integration  becomes  cumbersome  and  one  of  the  steps  in  the 
derivation  is  fallacious.  In  the  notation  of  Probert ’s  paper 
the  variables  are  t  and  d  the  diameter  at  time  t.  Drops  of 
diameter  d  may  be  considered  to  result  4  from  larger  drops,  /xi  -A 
injected  at  a  time  t  secs,  earlier,  of  diameter  at  injection  (o +XV) 

volume  of  drops  of  size  d  — >  in 

steady  str  ekis,  ^-4,  mV 

% 


the 


■V|8d  =  d  |  -*n 


e  W 


$(d  +At)  Vd  t 


t-=o 

Probert  puts  £d  in  place 
alternative  variables  d 


and  dx  ==  [a  - 


■xt) 


'h.  t  confusing  the 

This  expression 


can  be  integrated  over  values  of  d  from  zero  to  infinity  to  give 
the  same  result  as  (22)  but  the  evaluation  of  the  double  integral 
is  somewhat  involved. 
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/APPENDIX  B. 


APPENDIX  3. 


Galculntion  of  the  rate  of  combustion  of  mono  - 

propellant  droplets. 


Strictly,  to  evaluate  the  rate  of  combustion ,  vuc  need 
the  condition  for  which  a  solution  of  the  partial  different  ini 
equations  of  heat  flow  find  mass  flow  exists  under  the 
appropriate  boundary  conditions.  In  the  steady  state,  these 
equations  may  be  formulated  as  :- 

Heat  flow  balance;  div  (  \  grad  T) 

-  -a )  ■  A ’A  ~  c 

"X  being  the  thermal  conductivity 

the  temperature 

the  moss  flow  per  unit  time 

the  heat  content  per  unit  mass 

the  rote  of  reaction  per  unit  volume  of 
reactant  A  etc. 

the  heat  of  reaction  per  unit  mass  of 
reactant  A  etc. 

Mass  flow  balance:  div  (Da  grade  )  -  div  (Mn )  -+  W*  “O 

-  v/A  '  (2) 

being  the  diffusion  constant  of  reactant  A,  2  "being 
the  mass  fraction  of  A, 

E»  - 1 

For  multi -component  reactions,  the  mass  flow  equation  can 
be  repeated  for  each  component. 

Finally  conservation  of  mass:  div  (M)  =0  (3) 

and  wA  =  f  (a,  .  p,  T) ,  etc.  (4) 

the  form  of  the  function  depending  upon  the  kinetics  of  the 
reaction;  and  the  appropriate  boundary  conditions  for  T  and 
etc.,  and  their  derivatives.  Given  these,  and  if  such  n 
stationary  state  is  possible,  n  relation  exists  between  the  rate 
of  burning  and  the  other  parameters.  It  has  been  possible  to 
express  this  ns  an  algebraic  function  only  for  certain  cases,  few 
in  number,  and  then  only  with  the  help  of  drastic  assumptions. 

For  a  review  of  these,  reference  should  be  made  to  Corner  (39) 
and  Mnrkstein  and  Polanyi  (4r),  A  term  for  the  flow  of  heat 
by  radiation  might  be  included  in  equation  (l).  In  most  systems 
the  heat  transferred  by  radiation  is  small  compared  with  that  due 
to  conduction  and  the  term  is  usually  omitted.  In  combustion 
systems  where  the  velocities  of  flow  are  low  and  the  pressure 
variations  throughout  the  system  negligible,  the  translational 
kinetic  energy  can  be  neglected  and  the  equation  for  the  con¬ 
servation  of  momentum  omitted.  Boys  &  Comer  (ll)  have  calculated 
the  rate  of  burning  of  cordite  by  numerical  integration  of  the 
equations,  but  only  for  the  case  of  a  plane-  surface  of  propellant; 

/for 


T 

M 

H 

WA 

3-A 


-  div  (hH')+W*Qft=0 
(1) 
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for  the  problem  now  considered  the  spherical  symmetry  renders 
the  numerical  integration  much  more  laborious  because  the 
independent  variable  appears  explicitly  in  the  differential 
equations. 

The  difficulty  in  obtaining  algebraic  solutions  lies  in 
the  exponential  form  of  the  equation  connecting  reaction  rate 
and  temperature  which,  for  the  case  of  an  unimolecular  reaction 
is  : 

a  -e/rt 
w  s  apt)  e 

r  (5) 


Solutions  using  this  kinetic  form  have  been  obtained  for  a 
planar  unbounded  reaction  zone  by  Prank -Kamanet ski  (41)  but  the 
methods  used  do  not  appear  to  be  convenient  for  spherical  re¬ 
action  zones.  We  shall,  therefore,  employ  a  method  of  solution 
which  is  historically  the  oldest,  being  fundamentally  identical 
with  that  used  by  Mallard  and  le  Chatelier  (1883)  and  with  vari¬ 
ations  by  other  workers  up  to  the  last  war.  As  has  been  pointed 
out  *y  Corner  (loc.cit)  the  approximation  determines  the  variation 
of  rate  of  reaction  throughout  the  reaction  zone  and  solves  the 
heat  flow  equation  only  in  the  heating  zone  where  the  rate  is 
assumed  to  be  zero.  However,  if  we  are  content  with  obtaining 
the  rate  of  combustion  per  unit  surface  of  a  drop  of  diameter  d 
in  terms  of  the  rate  of  combustion  of  a  plane  surface,  the  reaction 
rate  can  be  eliminated  from  the  result.  The  parameters  necessary 
are  a  temperature  Tx  below  which  the  rate  of  reaction  is 
negligible  and  T0  the  temperature  of  the  liquid  surface.  These 
and  other  assumptions  are  discussed  on  p.  10,  11  &  IS.  If  the 
influence  of  diffusion  of  reactant  is  assumed  to  be  negligible 
equation  (S)  need  not  be  considered. 


For  a  spherical  drop  of  diameter  d  equation  (l)  becomes 
dT 


(6) 


T  &  w  being  the  temperature  and  reaction  rate  at  o  distance  r  from 
the  drop  centre. 

The  boundary  conditions  at  the  liquid  surface  are 


Md%(37)=ML  ,  T-T„ 

b  * 


(7) 


here  L  is  the  heat  of  vaporisation  in  calories/gram,  and  CD 
the  constant  pressure  specific  heat  in  calories/degree/gram- 
assumed  to  be  the  same  for  both  reactant  and  products. 

In  terms  of  the  dimensionless  variables  A  -  ^p(T-T0  +  Cp) 

L 

1  v  3* 


With  §=  MCp/^Xd 
I 

iqe 


(8) 


and  where  Cp  and  )\  are  given  mean  values  for  the  temoe return 
range  covered,  equation  (cT  may  be  written: 

OCfisTw 


d2e  , 


d  R 


4>  L  R4' 
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--0 


(9) 


/in 


UNCLASSIFIED 


in  the  region 


T0<T<T, 


W  •*  O 


Qnd  §* + Ur = 0 

do) 

with  boundary  conditions 

Po^'.Qo*1,  (sslr^ 

(n) 

the  solution  of  (10)  is 

io -m 

(12) 

0*  e 

In  the  reaction  zone  we  may  assume 
equation  (9)  becomes  - 

\d£i  -  ^Cpcl  .  w 
^dR  4>L  ^ 

cJG* 


$ 


to  be  constant,  then 
(13) 


To  eliminate  we  invoke  the  condition  that  0  must  be  o 

continuous  function  of  R  at  the  junction  of  the  heating  and 
reaction  zones,  then  nt  the  junction,  R  =  Rj_ ,  and  from  (12) 
and  (13) 


d0\  -  X0= 


but  from  (12)  R(  —  1 *" 
whence  (14)  becomes  ^  ?' / 1  _  =i  Q.Qp.j.  w, 

V  X  1  4>l  B\ 


(14) 


(15) 


If  we  assume  that  w;  >9;  are  independent  of  the  radius  of  the 
drop,  equation  14  is  a  relation  between  the  mass  r^tc  of  burning, 
the  drop  size  and  the  other  parameters  of  the  system.  Since 
r-'i  =s  Tfdafn  where  m  is  the  rate  of  burning  per  unit  surface, 
we  obtain  from  (15) 

fn^Cp  /.  .  2>  InOi^4  3 

>  V  J  LQ,  (16) 


as  d  tends  to  infinity  m  tends  to  mo  the  rate  of  burning  per 
unit  surface  for  a  plane  surface,  we  may  write 

.1/  2  A  Vi  0  A4-  -  m2 

^  V  "^riCpCT  /  (17J 


This  equation  for  rh  has  four  distinct  roots,  given  by  : 
2iri  =  ±Tn0  +  *Tp T  CfS  ) 


(18) 


toking  positive  values,  we  see  that  if  ^ 


If)  ^  Tre  ¥- 


Cpd 


(19) 
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/and 


UNCLASSIFIED 


and  if 


UNCLASSIFIED 


(sc.) 


In  view  of  the  assumptions  needed  for  the  derivation  of 
equation  (18)  it  cannot  be  expected  to  give  more  than  the 
order  of  the  variation  of  ih  with  diameter.  It  is  intended 
to  check  this  calculation  by  numerical  integration  of  equations 
(l)  and  (2)  when  sufficient  data  is  available  to  make  the 
task  profitable. 
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